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FOREWORD 


This  report  was  prepared  by  the  Computer  Applications  Section  of 
General  Dynamlcs/Electric  Boat  as  part  of  the  Submarine  Integrated 
Control  Program  (SUBIC)  of  the  Office  of  Naval  Research.  Electric 
Boat  Is  coordinator,  under  Contract  Nonr  2512(00),  of  this  program; 
Cdr.  P.  R.  Haselton,  Jr.,  USN,  is  Project  Officer  for  ONR;  Dr.  H.  E. 
Sheets  Is  Project  Coordinator  for  Electric  Boat;  and  Dr.  A.  J. 
van  Woerkom  Is  Chief  Computer  Scientist. 

The  program  Is  divided  into  several  parts:  ship  control,  weapon  and 
tactical  control,  engineering  control,  communications,  environmental 
control,  and  command  control.  This  report  is  one  of  a  series  dealing 
with  tactical  control. 


ABSTRACT 

(Unclassified) 

This  report  consists  of  a  mathematical  model  of  a  simple,  statistical 
fire  control  scheme.  Itie  equations  are  developed  for  the  determina¬ 
tion  of  relative  target  motion  parameters,  a  complete  bearings-only 
solution  and  complete  solutions  based  upon  hypothesized  Inputs  of 
target  speed,  course,  or  range.  The  statistical  properties  of  some 
of  the  results  are  analyzed  and  compared  to  Mark  113  results  and  Ideal¬ 
ized  manual  plots.  A  method  of  detecting  target  zigs  is  described 
and  statistically  evaluated.  Methods  of  improving  some  of  the  present 
Mark  113  fire  control  system  computational  schemes  are  presented. 
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SUMMARY 


Purposes 

Automatic  Statistical  Processing  (ASP)  is  a  fire  control  computational 
scheme  that  has  been  developed  to  provide  all  available  tactical  data 
rapidly  and  accurately  and,  ultimately,  to  provide  a  fire  control 
solution  for  a  target  employing  evasive  tactics. 

Techniques 

1 .  Bearing-time  C?urve  and  Relative  Motion  Analysis 

ASP's  basic  technique  consists  of  statistically  smoothing  bearings  to 
a  mathematical  approximation  to  the  bearing -time  curve.  By  using  a 
digital  computer  for  rigorous  statistical  processing,  the  best  possible 
curve  (least-square  fit)  is  produced.  The  smoothed  curve  used  in  man¬ 
ual  plotting  is  obtained  by  eye  and,  in  general,  only  approximates  the 
least -square  fit. 

The  equation  of  the  least-square  bearing-time  curve  yields  the  bearing 

•  •  • 

(B),  bearing-rate  (b)  and  bearing  acceleration  (B)  as  functions  of  time. 
These  parameters  contain  all  the  essential  Information  contained  in  the 
raw  bearing-time  data.  The  following  relative  motion  parameters  can  be 
deduced  from  this  Information: 

1.  relative  angle -on-t he -bow  (a) 

2.  relative  course  (O) 

3.  ratio  of  relative  speed  to  initial  range  (u/R  ) 

f  2 ) 

4.  minimum  target  speed 

In  order  to  obtain  a  meaningful  approximation  to  the  bearing-time 
ctirve  and  the  resulting  parameters,  it  is  necessary  for  own-ship  to 
travel  with  vinlform,  rectilinear  motion.  This  is  true  for  both  ASP 
and  manual  plotting  techniques. 
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ASP  also  has  the  capacity  to  produce  hypothetical  solutions  on  any 
single  leg  of  ovm-shlp  track  by  Introducing  values  of  target  speed, 
course  or  range  Into  the  computer.  The  values  Introduced  may  be 
estimated,  hypothesized  or  based  on  physical  limitations.  By  dis¬ 
playing  the  resulting  "solutions",  qualitative  Information  about 
possible  tactical  situations  Is  provided. 

2.  Zlg  Detection 

The  Initial  step  toward  ultimate  solution  of  the  passive  maneuvering 
target  problem  Is  the  detection  of  target  zlgs.  ASP  Includes  a 
technique  whereby  a  computer  calculates  the  probability  that  a  target 
has  zlgged.  The  principle  used  Is  fundamentally  the  same  as  that 
used  In  the  Barnard  plot.  A  least -square  bearing -time  curve  Is  ex¬ 
trapolated  to  obtain  an  estimate  of  the  expected  future  bearings. 

The  occurrence  of  a  target  zlg  is  then  indicated  by  systematic  de¬ 
viations  of  the  rneasxired  bearings  from  the  predicted  bearings.  Since 
deviations  will  result  from  bearing  noise  even  when  a  target  zlg  has 
not  occurred,  only  the  probability  of  a  target  zlg  can  be  determined. 
Probabilities  In  excess  of  9^  are,  hov/ever,  usually  Interpreted  as 
certainty.  Comparison  of  the  expected  bearings  to  the  measured 
bearings  Is  continued  for  approximately  two  minutes  whence  the  bearing 
data  obtained  In  the  tv.'o  minute  Interval  are  automatically  combined 
with  the  bearing  history  to  form  a  new  least -square  fit.  The  next 
two  minute  Interval  Is  then  scanned  In  the  same  manner.  Wien  a  zlg 
Is  detected,  the  history  prior  to  detection  Is  disregarded  (zlg  de¬ 
tection  reset)  and  a  new  history  Is  collected  on  the  next  leg  of 
target  track  for  the  purpose  of  detecting  possible  subsequent  target 
zlgs.  Since  the  zlg  detector  Is  also  sensitive  to  own-ship  zlgs.  It 
must  be  reset  when  own-ship  maneuvers. 

Two  types  of  bearlng-tlme  curve-fitting  have  been  Investigated  for 
purposes  of  zlg  detection;  a  curved  line  (quadratic  scheme)  and  a 
straight  line  (linear  scheme).  The  extrapolation  and  deviation 
caused  by  a  target  zlg  are  Illustrated  In  Fig.  C-1-1  on  page  116  for 
the  quadratic  scheme  and  In  Pig.  C-5-2  on  page  126  for  the  linear 
scheme. 
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3 .  Bearings -only  Complete  Solution 

When  the  zlg  detector  Indicates  that  the  target  Is  straight  running,  a 
complete  bearlngs-only  target  localization  solution  Is  possible.  To 
obtain  the  solution,  own-ship  travels  on  a  straight  leg  and  collects, 
bearlng-tlme  data,  then  changes  course  and/or  speed  and  collects  a 
second  set  of  bearlng-tlme  data.  By  combining  the  Information  contained 
In  the  resulting  two  statistically-smoothed  bearlng-tlme  curves,  the 
range,  course  and  speed  of  the  target  can  be  obtained. 

Results 

1.  Statistical  Error  Analysis 

Since  all  bearing  data  have  random  deviations,  all  parameters  derived 
from  bearlng-tlme  data  are  subject  to  statistical  uncertainty.  The 
curves  In  sections  3>  5t  and  6  of  Appendix  B  illustrate  the  \in- 
certaintles  in  the  relative  motion  parameters.  The  uncertainties  in 
bearing,  bearing-rate  and  bearing  acceleration  decreases  as  tracking¬ 
time  Increases,  The  uncertainties  in  the  other  relative  motion  para¬ 
meters  decrease  with  Increasing  bearing-rate  and  time  and  with  de¬ 
creasing  range.  Although  no  formal  error  analysis  has  been  made  for 
the  hypothetical  solutions,  their  accuracies  are  not  expected  to  be 
very  high.  However,  this  type  of  solution  can  be  used  to  furnish 
early  Information  about  possible  tactical  situations. 

The  percentage  error  in  bearlngs-only  range  decreases  as  tracking  times 
and  the  absolute  difference  In  the  bearing -rates  on  the  two  own-ship 
legs  Increase.  Statistical  evaluation  of  the  bearlngs-only  solution 
has  also  indicated  the  possibility  that  bearing-rate  and  its  uncertainty 
can  be  used  to  determine  the  optimum  time  for  own-ship  zlg.  The  de¬ 
tails  for  establishing  an  own-ship  zlg  time  criterion  are  yet  to  be 
determined. 

The  linear  zlg  detection  scheme  has  been  found  to  be  more  effective  for 
long  range  targets.  Figs.  C-6-1  and  C-6-2  and  the  accompanying  tables 
contain  quantatlve  comparisons  of  the  two  schemes. 


vil 


CONFIDENTIAL 


CONFIDENTIAL 


2.  Comparison  of  ASP  with  Manual  Plots 

To  quantatlvely  compare  ASP  with  manual  plots,  it  Is  necessary  to 
assume  that  manual  plots  are  rigorously  least-square  fitted.  Compari¬ 
son  can  then  be  made  directly  because  ASP's  mathematics  are  implicit 
in  manual  plotting  techniques. 

The  ASP  scheme  Is  superior  because  of  the  high  data  assimilation  rate 
possible  using  digital  computer  processing.  Parameters  are  obtained 
more  accurately,  more  rapidly  and  at  longer  ranges.  The  curves  of 
Section  6  of  Appendix  B  (pages  82  through  106)  compare  the  accuracies 
achieved  by  the  two  techniques  for  various  geometries.  A  complete 
discussion  of  the  assumptions  and  results  of  the  comparison  is  con¬ 
tained  In  Sections  4. A  and  4.B  of  Part  II  of  this  report  (pages  29 
through  31)* 

3 .  Comparison  of  ASP  with  Mark  113 

The  Mark  113  system  Includes  relative  motion  analysis  and  a  complete 
bearlngs-only  solution  (Churn). 

Relative  motion  analysis  In  the  present  Mark  113  Is  similar  to  that  In 
ASP  In  that  bearing-time  data  are  used  as  Input  and  own-ship  Is  con¬ 
strained  to  uniform,  rectilinear  motion  to  obtain  a  meaningful  relative 
motion  analysis.  However,  the  only  parameters  provided  by  the  Mark 
113  relative  motion  analysis  are  the  Initial  angle-on-the-bow  and  the 
ratio  of  relative  speed  to  Initial  range.  Since  no  statistical  error 
analysis  has  been  made  for  the  Mark  113  relative  motion  parameters, 
accuracies  for  the  two  techniques  cannot  be  compared. 

In  the  Mark  113  system,  a  hypothetical  solution  can  be  obtained  for  an 
Input  of  target  speed  only.  In  the  event  that  a  dual  .j^^lutlon  exists 
(see  Fig.  A-4-2  c  page  51 )>  only  the  longer  range  solution  Is  pro¬ 
vided,  whereas  ASP  furnishes  both  solutions. 

The  Mark  113  bearlngs-only  solution  uses  bearing-time  data  and  own- 
ship  motion  parameters  for  at  least  two  own-ship  legs  to  provide  tar¬ 
get  range,  course  and  speed.  Although  Churn  theoretically  allows 
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own-ship  freedom  of  motion  (with  resulting  mathematical  complexity), 
in  reality,  an  optimum  tactic  (two  straight  legs  differing  In  course 
and/or  speed)  Is  necessary  to  produce  adequate  solutions  at  ShhiiOC 
ranges.  Also,  own-ship  straight  legs  are  required  in  order  to  detect 
target  zlgs.  It  is  therefore  apparent  that  the  freedom  of  own-ship 
motion  provided  by  the  Mark  113  system  Is  not  a  requirement  In  connec. 
tlon  with  SUBROC  fire  control  solutions. 

4.  Bearing  Pre -Smoothing 

In  addition  to  being  complex,  the  machematlcs  In  Churn  Include  an  In¬ 
herent  bias.  It  has  been  shown  (reference  l)  that  this  bias  can  be 
reduced  to  insignificance  by  pre-smoothing  (averaging)  groups  of  bear 
mgs.  At  present  the  averaging  criterion  is  based  on  bearing-rate. 
The  ASP  study  Indicates  that  a  better  criterion,  based  on  bearing 
acceleration  and  its  uncertainty,  can  be  developed. 
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INTRODUCTION 


Automatic  Statistical  Processing  (ASP)  Is  a  continuation  of  earlier 
SUBIC  studies  in  the  fire  control  area.  The  objective  set  forth  for 
ASP  (in  addition  to  the  more  general  SUBIC  objectives)  Is  to  provide 
all  useful  tactical  lnfori..atlon  which  can  be  derived  from  available 
data  as  rapidly  and  accurately  as  possible. 

It  Is  felt  that  a  significant  portion  of  available  tactical  Informa¬ 
tion  Is  not  now  displayed  In  present  automatic  computing  systems. 

For  example.  In  the  existing  Mark  113  fire  control  system,  the  Initial 
relative  angle -on- the -bow  and  ratio  of  the  relative  target  speed  to 
Initial  range  are  the  only  quantities  derived  from  the  bearings  on  the 
first  leg  of  own-ship  maneuver.  Range,  course,  and  speed  are  avail¬ 
able  either  upon  entering  a  target  speed  estimate  or  after  an  own-ship 
zlg.  Other  tactical  Information  such  as  target  zlg  indications, 
smoothed  bearings  and  bearing -rate ,  a  quantity  often  used  to  fire 
acoustic  torpedoes.  Is  obtained  from  manual  plots  which  are  severely 
limited  in  accuracy  at  the  longer  ranges . 

Ihe  technique  used  to  meet  the  ASP  objective  consists  of  smoothing. 

In  the  statistical  sense,  bearing  data  to  a  mathematical  approxima¬ 
tion  to  the  bearing -time  curve.  This  technique  provides  relative  mo¬ 
tion  parameters  during  the  first  leg  of  own-ship  maneuver,  results  In 
simplified  solution  equations,  and  makes  zlg  detection  possible.  In 
order  to  achieve  these  ends,  however.  It  Is  necessary  to  constrain 
own-ship  to  uniform  rectilinear  motion.  Relative  motion  parameters 
can  be  obtained  on  any  own-ship  uniform  rectilinear  motion  leg.  Such 
legs  will  be  referred  to  as  single  legs  throughout  this  report;  and 
relative  motion  analysis,  or  simply  single  leg  analysis  will  be  used 
to  designate  calculations  leading  to  relative  motion  parameters. 

Ilie  relative  motion  parameters  which  can  be  obtained  from  any  single 
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leg  analysis  are: 

1.  relative  angle -on -the -bow 

2.  relative  course 

3.  statistically  smoothed  bearings 

4.  bearing -rate 

5.  change  In  bearing-rate 

6.  ratio  of  relative  speed  to  Initial  range 

7.  mlnlmxim  target  speed 

niese  quantities  can  be  supplied  continuously  for  the  entire  tracking 
period.  It  Is  also  possible  for  the  operator  to  introduce  Into  the 
computer  estimates  or  hypothetical  values  of  target  range,  course  or 
speed  to  obtain  possible  target  localization  solutions.  Ihe  possible 
solutions  can  be  geographically  represented  on  a  display  similar  to  a 
strip  plot  for  a  qualitative  "picture"  of  the  possible  tactical 
situations . 

When  confronted  with  a  zigglng  target,  as  indicated  by  the  zlg  detect¬ 
or,  It  may  be  possible  to  obtain  a  range  estimate  by  obtaining  para¬ 
meters  for  each  leg  of  the  target  track  and,  from  these,  selecting 
portions  of  the  possible  target  tracks  where  zigs  will  not  severely 
degrade  range  solutions.  In  the  event  that  the  target  Is  slnuatlng. 

It  may  be  possible  to  approximate  the  track  by  a  series  of  straight 
legs  and  use  the  same  technique  as  for  a  zigglng  target. 

Thus,  a  leg -by -leg  analysis,  coupled  with  the  ability  of  a  trained  op¬ 
erator,  might  possibly  be  used  to  partially  overcome  the  difficulties 
presented  by  a  zigglng  target  so  that  a  rough  estimate  of  target  range 
can  be  passively  obtained.  This  range  could  be  used  to  Indicate  whether 
the  target  Is  within  weapon  range  and,  if  so,  whloh  weapon  to  use. 
Single  ping  techniques  could  then  be  used  to  obtain  a  more  accurate 
range.  If  necessary.  Hius,  the  possibility  of  pinging  a  target  beyond 
weapon  ranges  could  be  reduced. 
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In  addition  to  the  objectives  and  capabilities  mentioned  above,  the 
ASP  fire  control  system  can  also  be  used  to  obtain  a  complete  bearlngs- 
only  solution.  To  obtain  this  solution,  own-ship  maneuver  must  consist 
of  two  straight  legs  differing  In  course  and/or  speed.  Comparison  of 
the  bearlngs-only  solution  to  the  hypothetical  solutions  obtained  on 
first  leg  of  own-ship  maneuver  can  serve  to  Increase  operator  confi¬ 
dence  In  computer  outputs . 

This  report  Is  preliminary  and  represents  the  results  of  the  study  to 
date.  Conclusions  based  on  this  work  and  an  outline  of  the  work  con¬ 
sidered  for  the  future  are  contained  In  Part  IV. 
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I 

TARGET  LOCALIZATION  SOLUTIONS 


1 .  Geometric  Relations,  Definitions,  and  Symbols 

•nie  expression  target  localization  Is  used  here  to  designate  the  data 
processing  which  provides  Information  about  the  time -dependent  posi¬ 
tion  of  the  target  relative  to  own-ship.  Complete  knowledge  of  the 
target's  relative  position  as  a  function  of  time  (subject  to  statisti¬ 
cal  uncertainties)  will  be  called  a  complete  solution.  TSius,  target 
range  and  bearing  as  functions  of  time;  or  target  range,  course  and 
speed,  and  own-ship  motion  constitute  a  complete  solution. 

T^ie  geometry  for  an  earth-fixed  coordinate  system  Is  shown  In  Pig. 
I-lrl  where  both  own-ship  and  target  are  traveling  with  uniform  mo¬ 
tion*.  Pig.  1-1-2  depicts  the  same  situation,  but  with  the  coordinate 
system  origin  fixed  to  ov/n-ship.  The  velocity  vector  relations  are 
illustrated  In  Pig.  1-1-3. 

The  symbols  used  In  the  aforementioned  figures  are  also  used  In  the 
mathematical  development.  Vector  quantities  are  designated  by  horiz¬ 
ontal  bars  above  the  symbols;  their  magnitudes  are  designated  by  the 
symbols  without  the  bars.  Superscripts  (l)  and  (2)  refer  to  own-ship 
and  target  quantities,  respectively.  Inputs  to  the  target  localiza¬ 
tion  computer  other  than  bearing  data  and  own-ship  parameters  are 
designated  by  an  e  subscript. 

Equations  appearing  In  Sections  I,  II,  and  III  of  this  report  are  de¬ 
veloped  In  the  appendices  and  bear  the  same  numbers.  The  useful  equa¬ 
tions  for  target  localization  solutions  are  enclosed  In  rectangles. 


*Thls  will  be  assumed  throughout  unless  otherwise  Indicated. 
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2.  BearIng-TIme  Relation  and  Approximation 

Topical  bearing -time  relations  are  qualitatively  represented  by  Pig. 
1-2-1  for  U/R^  both  small  and  large.  Ihe  point  of  maximum  bearing- 
rate  Is  the  closest  point  of  approach  (CPA)  of  the  target. 

It  Is  shown  In  Section  1  of  Appendix  A  that  limited  portions  of  the 
bearlng-tlme  relation  can  be  approximated  by 


In  which 


and 


B  =  a+b(t-tQ)+c(t-tQ)^ 


b  =»  (U/i^o)  3ln  a^, 

2 

c  =  sin  cos  a  . 

'  o  o  o 


(A-1-2) 


(A-1-3) 

(A-1-4) 

(A-1-5) 


The  approximation  to  the  bearing  time  curve  (equation  A-1-2)  Is  more 
accurate  at  long  ranges  than  at  short  ranges . 

A  least-square  processing  of  bearings  versus  time  using  (A-1-2)  yields 
the  matrix  equation 


In  which 


(A-1-9) 


(A-1-10) 


I  ■ 
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tr 

V  = 

(A-l-ll) 


(A-1-12) 


and  is  the  bearing  obtained  at  time  t^ . 

Ihe  coefficients  a,  b,  and  c  as  obtained  from  equations  (A-1-9,  10) 
are  used  in  the  subsequent  section  to  obtain  relative  motion  pareuneters. 

3.  Relative  Motion  Parameters 

The  relative  motion  parameters  are  those  quantities  which  are  obtain¬ 
able  on  a  single  leg  of  own-ship  maneuver. 

Itie  equations  of  this  section  are  developed  in  Sections  2  and  3  of 
Appendix  A . 

Differentiating  equation  (A-1-2)  with  respect  to  time  gives 
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B  =  b  +  2c(t-t  ) 
o' 

which  when  differentiated  again  gives 


•  • 

B  =  2c 


(A-2-1) 


(A-2-2) 


B  and  B  are  bearing-rate  and  the  change  in  bearing-rate,  respectively. 

The  change  in  bearing-rate  Indicates  whether  the  target  is  opening  or 
closing. 


The  relative  angle -on-the -bow  is  obtained  from  the  relation 


Oq  =  tan-^(bVc) 


for  the  initial  value  and  from 


a(t)  =  a^+B(t)-a 


(A-3-1) 


(A-3-2) 


for  any  subsequent  time  (t). 


Relative  course  can  then  be  found  from  the  relation 


0  “  TT  -  +  a 


(A-3-3) 


The  ratio  of  relative  speed  to  initial  range  can  also  be  obtained 
from  bearing  data  alone.  The  relation  is  expressed  by 
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or 


U/Rq  =  b/sin 


(A-3-4) 


r  2  4 1 

1/2 

c  +b 

(A-3-5) 


4.  Requirements  for  a  Complete  Solution 

Bearing  data  alone  Is  Insufficient  to  obtain  a  complete  solution  when 
own-ahlp  travels  with  uniform  motion.  It  is  possible,  however,  to  ob¬ 
tain  a  complete  solution  when  either  range,  course,  or  speed  of  tar¬ 
get  Is  specified  (l.e.,  obtained  externally).  Occasionally,  estlniates 
of  one  or  more  of  these  parameters  are  available  and  can  be  used  to 
obtain  estimates  of  the  remaining  parameters.  Alternatively,  values 
can  be  hypothesized  In  order  to  obtain  various  possibilities  for  the 
remaining  parameters.  Thus,  for  example,  a  submariner  might  ask  for 
possible  ranges  and  courses  under  the  hypotheses  that  the  target  speed 
Is  ten,  twelve,  and  fifteen  knots.  Hypothesized  parameters  are  per¬ 
haps  more  prevalent  than  estimated  parameters  In  the  manual  and  semi¬ 
automatic  plotting  systems. 

5.  Solution  Using  Speed  Input 

In  this  and  the  following  sections.  Inputs  to  the  target  localization 
computer  other  than  bearing  data  and  own-ship  parameters  are  desig¬ 
nated  by  an  e  subscript. 

•Rie  equations  of  this  section  provide  a  complete  solution  from  bearing 

(2) 

data  and  a  given  target  speed  ' ) .  The  development  Is  In  Section  4 

of  Appendix  A. 
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(A-4-2) 


Equation  (A-4-2)  indicates  the  possibility  of  two  solutions.  This 
dual  solution  possibility  Is  a  property  of  the  fire  control  geometry 
and  Is  Independent  of  the  particular  mathematical  descriptions  used 
(see  Section  4  of  Appendix  A).  It  is  also  possible  to  enter  a  speed 

which  Is  too  small  to  provide  a  solution.  ’This  latter  situation 
leads  to  a  minimum  target  speed  (subject  to  statistical  uncertainties). 

By  examining  the  angle  ()2^^  -a+a^)  and  the  various  solution 

possibilities  may  be  determined.  Fig.  1-5-1  summarizes  the  different 
possibilities . 

When  more  than  one  solution  exists,  two  values  of  y  exist.  One  Is 
given  by  the  principal  value  of  equation  (A-4-2)  and  the  other  by  Its 
complement .  Thus 


II 

(A-4-4) 

(A-4-5) 
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In  the  event  that  tv;o  y  angles  are  possible,  equation  (A-4-1)  and 
those  which  follow  v/111  have  to  be  duplicated  to  provide  the  tv^o  pos¬ 
sible  solutions.  Thus, 


=  TT-Y.+a-a^  (J  =  1  or  1,2) 


(A-4-7) 


U, 


,  (A-4-8) 


,  and 


(A-4-9) 

and 

(A-4-10) 


in  which  R(t)  Is  the  range  at  any  time  t. 

In  many  cases,  one  of  the  two  solutions  v/lll  be  physically  absurd  al¬ 
though  It  Is  mathematically  possible. 


6.  Solution  Using  Range  Input 

The  equations  of  this  section  are  derived  In  Section  5  of  Appendix  A. 
Let  RgCtg)  be  a  range  Input  at  time  t^,  then 


Vj^^^^(te)  =  V^^^sin 

b+2c(te-to)] 

(A-5-6) 
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■5 


and 


cos 


0 


(1)_ 


B(tJ 


Re(te)c 


in  which  V  )  and  V  are  the  cross -llne-of -sight 

X  e  ye 

llne-of -sight  speeds  at  t^ ,  respectively,  and 


n(t^)  =  a4b(tg-t^)^c(tg-t^)^ 


Target  speed  Is  then  obtained  by 


(te)^Vy 


(2)' 


1 1/2 


("e) 


and  the  target  course  by 


=  P(tg)4tan 


(s) 

^7^ 


Range  at  any  time  other  than  t^  is  then 


R(t) 


R^(te)  [b-+c(t^-t^)] 
b+c(t-tQ) 


IT 


(A-5-7) 

and  along- 

■(A-5-2) 

(A-5-8) 

(A-5-9) 

(A-5-10) 
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7.  Solution  Using  Course  Input  (See  Section  6  of  Appendix  A  for  development) 

f  2  ^ 

For  an  estimate  of  tarf^et  course  (0^'  '),  target  speed  is  obtained 
from 


in  which 


(A-6-1) 


T 

I 

I 

I 

I 

I 

I 

I 


Target  range  can  then  he  found  fr-oin 


8.  Bearings-Only  Solution  (See  section  7  of  Appendix  A  for  development) 

For  a  bearings-only  solution,  own-ship  travels  for  a  time  t'  with  uni¬ 
form  motion,  then  changes  course  and/or  speed  and  travels  with  uniform 
motion  until  the  end  of  the  tracking  time.  The  simplifying  assumption 
is  made  that  own-ship  changes  course  and/or  speed  instantaneously  at 
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t*.*  Bearing  data  is  collected  and  processed  on  each  of  these  two 
"legs"  of  own-ship  track.  The  subscripts  1  and  2  are  used  to  Indicate 
first  and  second  leg  parameters,  respectively.  The  x  and  y  subscripts 
are  used  to  indicate  the  cross -llne-of -sight  and  along-llne-of -sight 
speeds,  respectively. 

At  the  time  of  own-ship  zig  (f),  it  can  be  seen  that  two  values  of 
own-ship  cross-line-of -sight  speed  and  two  values  of  bearing-rate  are 
obtainable,  one  from  each  of  the  legs  of  own -ship  track.  Thus,  for 
leg  1, 


Vix^^^f)  =  V^^^^sin[jZf^(^^-B(t')] 


(A-7-4) 


and 


B^(t')  =  b^+2c^(t'-t^) 


(A-7-7) 


For  leg  2, 


(A-7-5) 


and 


(A-7-8) 


♦See  Section  2,  Part  IV. 
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The  bearing  B(t')  is  common  to  both  legs  and  is  calculated  using  the 
afjproprlate  time  weighing  factors.  Thus, 


t-t* 

t 


^2 


(A-7-6) 


llie  range  at  t'  Is  then  calculated  from 


R(f) 


B2(fc')-B^(t') 


(A-7-3) 


and  for  any  subsequent  time  (t). 


R(t) 


R(f)b2 

bg+CgCt-t ' ) 


(A-7-9) 


Target  speed  and  course  are  constant  by  assumption  and  can,  therefore, 
be  calculated  using  Information  obtained  on  both  legs  of  own-ship  man¬ 
euver.  Thus,  for  leg  k  (k  =»  l,  2), 


,  (A-7-10) 


and 


(A-7-11) 
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in  which  B(t')  Is  found  by  equation  (A-7-6), 


R(t')CT 


and 


Using  the  above  equations. 


O 

* 

1 1 

and 


=  B( t ' )+tan”^ 


Introducing  the  appropriate  weighing  factors. 


r(2) 


W  V 


(2) 


WgVg 


(2) 


and 


(A-7-12) 


(A-7-13) 


(A-7-14) 


(A-7-15) 


(A-7-16) 


(A-7-17) 
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11 

STATISTICAL  PROPERTIES 


1.  Variances  and  Covariences  of  Least-Square  Coefficients  (see  Section  i  of 

Appendix  B  for  development) 

The  variances,  covariances,  and  correlation  coefficients  of  the  coef¬ 
ficients  a,  b,  and  c  are  ,  Ivon  oy 


(B-1-2) 


(B-1-3) 


(B-1-4) 
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where  the  c^^'s  are  the  elements  of  the  inverse  of  matrix  A,  v;hlch  Is 
obtained  from  the  coefficients  of  a,  b,  and  c  of  equations  (A-1-7)  or 


A  =  2(t^-t^) 


i(\-to) 


The  elements  of  A  consist  of  somev/hat  cumbersome  summations.  Tbere- 
fore,  to  facilitate  calculations  for  analysis  it  is  assumed  that  the 
bearings  are  obtained  at  equal  time  Increments  (T).  The  variances 
then  become 


2  9  2 

°a  ^  n°  V,  > 


o  ^  o 


{B-l-12) 


o  2  ^  1^  o  2 

c  TTs  E  • 


In  which  Og  is  the  standard  deviation  of  sonar  bearings  and  n  Is  the 
ntimber  of  bearings  sampled. 


The  covariances  are 


T  ,  O  O,  = 

ao  a  b 


36  „  2 
7,^  B  • 


fv.  0^0  = 

be  b  c 


180  „  2 


(B-l-13) 
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TOO 

ac  a  c 


T^n-" 


and  the  correlation  coefficients  are 


^ab  =  =  -  0.866, 

"'be  “  =  -  0.97,  (B-1-14) 

^ac  ^  VV3  =  0.749. 

Pigs.  B-1-1,  B-1-2,  and  B-1-3  of  Section  1,  Appendix  B,  are  plots  of 
the  standard  deviations  of  the  least  square  coefficients  as  functions 
of  time.  A  bearing  sampling  interval  of  two  seconds  and  a  sonar  bear¬ 
ing  standard  deviation  of  0.2°  are  used. 

2.  Variance  in  Relative  Motion  Parameters  (See  Sections  3  and  4  of  Appendix 
B  for  development) 

The  standard  deviation  in  bearing  rate  is  given  by 


(B-3-2) 


The  standard  deviation  in  the  change  in  bearing  rate  is  given  by 


(B-3-3) 


The  standard  deviation  in  a  is 

o 


c 


(B-4-2) 
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o  =  o 
a  a 


(B-4-3) 


Also,  the  standard  deviation  In  0  Is 


»c  =  =% 


(B-4-7) 


The  general  expression  for  standard  deviation  In  the  parameter  U/R^ 
Is  given  by 


'U/So 


,  2  2 
Sin  a  -cos 
_ o _ 

sin  a. 


b  -  (U/Rq)  cot  a  Oj  .  (B-5-1) 


For  a  S'  Tr/2,  equation  (B-5-1)  reduces  to 


o  ,  0 

U/1^o 


(B-5-2) 


For  »  0, 


(B-5-3) 


but  In  this  case  bearing  rate  Is  zero  and  the  only  Information  which 
can  be  derived  from  the  data  Is, 

B  =  constant 


B  =  0. 
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The  figures  in  Sections  3^  and  5  of  Appendix  B  illustrate  the  var¬ 
iation  of  the  standard  deviations  of  the  relative  motion  parameters 
'with  tracking  time. 

3.  Bearings-Oniy  Range  Variance 

Tiie  equations  of  this  section  are  developed  in  Section  7  of  Appendix  B. 


The  variance  of  the  range  at  current  time  (t)  is  obtained  from 


°Rlt 

I  = 

^  \** 

lA 

T 

8  V~3  T^/^o 


B 


I  A  B| 


(B-7-3) 


and  the  variance  of  the  I’ange  at  zir  time  (f)  is  obtained  from 


iA‘  8^ 

b| 


(B-7-6) 


in  which  T^  is  the  total  tracking  time  and  the  1  and  2  subscripts  re¬ 
fer  to  first  and  second  leg  quantities,  respectively.  The  above  two 
equations  show  that  the  range  uncertainty  decreases  with  the  total 
tracking  time  and,  for  a  given  total  time,  is  minimized  by  maximizing 
the  change  in  bearing-rate  (AB) . 


The  change  in  bearing-rate  can  be  very  nearly  maximized  by  running 
own-ship  perpendicular  to  the  initial  bearing  line  for  a  period  of 
time,  and  running  in  the  opposite  direction  for  the  remainder  of  the 
tracking  time.  The  terms  v/lthin  the  brackets  indicate  how  the  ratios 
of  the  times  on  each  leg  to  the  total  time  affect  the  range  uncertainty. 
The  bracketed  part  of  equation  (B-7-6)  is  a  minimum  vihen  T^  =  ^2 .  The 
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mlnimxam  for  the  bracketed  part  of  equation  (B-7-3)  Is,  however,  a 
function  of  both  the  times  and  the  bearing-rates  on  each  leg.  When 
the  bearing-rate  on  the  second  leg  approaches  zero,  the  uncertainty 
of  the  current  range  Increases  beyond  bound.  This  Is  also  a  property 
of  the  Mark  113  (Mode  2)  range  solution  as  time  on  second  leg  In¬ 
creases  beyond  bound  (Reference  1),  A  small  Bg  can  be  avoided  by 
running  own- ship  In  the  direction  of  the  bearing  drift  on  the  first 
leg  and  against  the  bearing  drift  on  the  second  leg.  If  for  some 
reason  this  Is  not  done  and  Bg  Is  small  compared  to  Bj^,  R(t'  )  will  be 
a  much  more  reliable  estimate  than  R(t),  The  Mark  113  system  com¬ 
putes  only  R(t),  (Reference  2), 

Ibe  minimum  for  the  bracketed  expression  of  equation  (B-7-3)  and  the 
comparable  expression  for  Mode  2  is  illustrated  in  Pig.  B-7-1  of  Ap¬ 
pendix  B.  The  X  and  y  variables  are  defined  by 


I 

I 

I 

I 

T 

I 

I 

I 

I 


The  ratio  Tg/T^  which  minimizes  the  braclccted  expression  of  equation  ^ 

(B-7-3)  is  designated  by  x-opt  and  represents  the  optimum  time  ratio 

for  a  given  y.  Positive  y  values  correspond  to  Bg  B^.  j 

Ibe  cusp  of  the  ASP  curve  shown  on  Fig.  B-6-1  does  not  extend  to  a  j 

•  I 

zero  value  for  x-opt  since  at  this  point  Tg  and  B^  would  also  be  zero, 
causing  the  bracketed  expression  in  equation  (B-7-6)  to  become  inde^  'i 

terminate.  This  expresses  the  fact  that  no  solution  can  be  obtained 
on  a  single  leg  (l.e.,  Tg  /  O).  .j 

The  minimum  value  of  the  bracketed  expression  for  various  values  of  y 
is  designated  by  z-opt.  Pig.  B-7-2  of  Section  7,  Appendix  B,  lllus-  | 

trates  z-opt  for  both  ASP  and  Mode  2.  In  general,  the  ASP  range 

;i 
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uncertainty  v;lll  be  about  1  Ji  times  la'r^or  than  the  Mode  2  uncertainty, 
but  when  approaches  zero,  the  range  error  will  be  the  same  for  both 
systems. 

The  term  optimum  tactic  Is  defined  as  the  ov:n-3hlp  maneuver  v/hlch  will 
result  In  the  smallest  range  uncertainty.  Optimum  tactic  designates 
both  the  direction  of  travel  and  the  length  of  time  on  each  leg.  Itie 
difference  In  the  optimum  tactics  for  ASP  and  Mode  2  lies  only  In  the 
length  of  time  on  each  leg.  For*  ASP,  the  time  on  the  first  leg  should. 
In  general,  be  slightly  more  than  half  the  total  tracking^ time  (T^  = 

.55  T^;  T2  "=  •  ^^5  Mode  2  (Reference  l)  the  time  spent  on 

the  first  leg  should  be  only  about  thirty  per  cent  of  the  total 
(T^  =  .3T^;  To  =  .TT^) . 


4.  Examples  for  Particular  Geometries  and  Comparisons  of  ASP  with  Manual  Plots 
A.  Relative  Motion  Parameters 


Itie  figures  of  Section  6,  Appendix  P,  show  relative  motion  parameter 
accuracies 


o 


a  ’ 


°"/«o 


o  • 

B 


and 


B 


f I.  both  ASP  and  manual  plotting  as  functions  of  time  for  the  partlcu- 
lai’  geometries  Illustrated.  Although  bearing -rate  Is  the  only  one  of 
these  parameters  presently  obtained  from  the  manual  plots,  the  mathe¬ 
matical  approach  for  obtaining  all  of  the  parameters  Is  Implicitly 
contained  In  plotting  techniques.  This  fact  is  used  to  make  a  compar¬ 
ison  of  man  versus  computer.  The  comparison  is  idealized  by  assvunlng 
that  the  manual  plots  represent  least -square  fits  of  the  data.  The 
essential  difference  between  the  man  and  the  computer  Is  taken  to  be 
the  data  assimilation  rate.  One  bearing  per  two  seconds  Is  assumed 
for  ASP  while  one  bearing  per  minute  Is  assiimed  for  the  manual  plots. 
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All  sonar  bearing  errors  are  considered  Independent  and  normally  dis¬ 
tributed  with  a  standard  deviation  of  0.2  degrees;  but,  since  physical 
plotting  of  data  Introduces  some  error,  the  total  standard  deviation 
of  bearing  data  for  the  manual  plots  Is  assiuned  to  be  0.3  degrees. 

For  these  sampling  Intervals  bearing-error  correlation  times  of  tv;o 
seconds  or  less  will  not  degrade  the  results  (Reference  1).  Likewise, 
a  constant  bearing  offset  will  not  degrade  the  results  although  the 
final  least-square  bearings  will  be  offset  by. the  same  amount. 

Itiese  graphical  comparisons  Illustrate  that 

1.  The  standard  deviations  In  the  relative  motion  parameters 
are  always  approximately  eight  times  less  for  ASP  than 
for  manual  plotting  when  the  same  geometry  and  tracking 
time  are  assumed. 

2.  To  achieve  equal  accuracy  with  manual  plotting  techniques 
longer  tracking  times  are  required. 

a.  For  the  required  time  using  manual  plots  Is  more 
than  four  times  as  great  as  that  required  for  ASP. 

b.  For  and  og  the  required  time  Is  about  two  and  one- 
thlrds  greater  than  that  required  for  ASP. 

c.  Since  is  geometry  dependent,  the  required  time 

varies  from  two  and  one -third  to  about  four  times  as 
great  as  that  required  for  ASP. 

B.  Manual  Plot  vs.  ASP  Range  Accuracy 

For  the  same  geometry  and  length  of  tracking  time,  the  standard  devia¬ 
tion  of  bearings -only  range  for  the  best  least-square  manual  plot  Is 
approximately  8.2  times  as  great  as  the  standard  deviation  obtained 


.1 

I 

] 

j 
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employing  the  ASP  system.*  This  difference  In  accuracy  Is  a  result 
of  the  ASP  system's  faster  bearing  saimpllng  rate  and  freedom  from 
errors  Incurred  In  physical  plotting. 

To  obtain  a  manual  plot  range  accuracy  comparable  to  that  obtained  by 
the  ASP  solution  the  necessary  tracking  time  Is  more  than  four  times 
greater  than  that  required  for  ASP. 

The  ratio  ° (see  figures  of  Section  6  of  Appendix  B)  gives  the  sig¬ 
nificance  of  bearing-rate  as  a  function  of  tracking  time.  Since  all 
bearings -only  solutions  are  dependent  on  the  change  In  bearing-rate 
resulting  from  an  own-ship  zlg,  the  range  uncertainty  will  be  a  func¬ 
tion  of  for  each  leg.  This  provides  a  basis  for  determining  the 

length  of  time  own-ship  should  track  on  a  single  leg,  the  details  of 
which  have  not  yet  been  determined.  No  such  criterion  Is  given  In  the 
present  Mark  113  system. 


I 


5.  Bearing  Pre-Snnoothing 


The  ratio  Indicates  the  significance  of  curvature  In  the  bearlng- 

tlme  curve.  It  has  been  shown  (Reference  1)  that  the  Inherent  bias  In 
the  Mark  113  bearings -only  solution  can  be  reduced  to  Insignificance 
by  averaging  ( pre -smoothing)  groups  of  bearings  and  using  the  result¬ 
ing  weighed  averages  as  single  bearings  for  the  least-square  process¬ 
ing.  No  Information  Is  lost  by  this  scheme  provided  the  bearlng-tlme 
curve  has  no  significant  curvature  within  the  averaging  Intervals. 

The  present  Mark  113  criterion  for  establishing  the  widths  of  the 
averaging  Intervals  (Reference  2)  is  based  on  the  bearing -rate .  This 


criterion  should  be  based  on  the  ratio. 


A  valid  criterion  Is  obtained  by  establishing  a  level  of  significance 
(K)  for  the  curvature.  Thus,  when 


*Uslng  the  assumptions  of  the  preceding  paragraph. 
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°B/6'  ^  ^  (K  a  constant) 

significant  curvature  Is  present.  The  first  qn  bearings  are  averaged 
(q  <  1  and  n  Is  the  total  number  of  bearings)  and  the  least-square 
process  Is  applied  to  the  remaining  bearings  and  the  Incoming  bearings 
\intll  Ogyg  >  K  when  the  entire  process  Is  repeated. 

Since  It  Is  necessary  to  have  at  least  four  bearings  for  a  solution, 
an  upper  limit  la  set  for  the  smoothing  time.  'Rius,  If 

B/fe 

when  the  time  Interval  Is  N  minutes,  the  bearings  In  the  N  minute 
Interval  should  be  averaged  and  a  new  Interval  started. 


I 
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III 

LONG  RANGE  BEARINGS-ONLY  ZIG  DETECTibN 


1 .  Detection  Scheme 

All  present  bearings-only  range  finding  methods  are  based  on  the  as¬ 
sumption  that  the  target  track  consists  of  uniform  rectilinear  motion. 

•nius,  the  usefulness  of  these  bearings-only  solutions  depend*  on  the 
ability  to  detect  target  zlgs. 

A  zlgging  target  Is  defined  here  as  one  which  does  not  conform  to  the  ' 

aforementioned  hypothesis.  Ibis  definition  includes  sinuating  targets.  ' 

( 

Hie  equations  of  this  section  represent  a  zlg-detection  technique 
v;hich  is  fundamentally  the  same  as  that  used  in  the  manual  plots.  In 
this  technique  bearing  data  are  collected  for  a  time  interval  and 
least-square  fitted  to  a  function  of  time.  Estimates  of  the  expected  ' 

future  bearings  are  obtained  from  the  least-square  fit  (bearing  his¬ 
tory).  The  occurrence  of  a  target  zlg  is  then  indicated  by  systema¬ 
tic  deviations  of  the  measured  bearings  from  the  predicted  bearings. 

When  systematic  deviations  are  not  discernible,  the  bearing  measure¬ 
ments  obtained  in  the  prediction  interval  are  combined  with  the  his¬ 
tory  bearings  to  form  a  new  history  from  which  a  new  least-square  fit 
is  obtained.  Bearings  are  then  predicted  for  the  next  interval  and 
the  entire  process  is  repeated.  Hie  entire  bearing -time  curve  is 
scanned  in  this  manner. 

2.  Equations 

Hie  tracking  time  is  divided  into  intervals  t^,  t^,  ...  tj,  *** 

such  that 


( t t **  At 
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for  (1  =  1,2,...).  Each  At  interval  contains  p  bearing  sampling 
points.  The  absolute  value  of  the  sum  of  the  deviations  of  the  pre¬ 
dicted  bearings  from  the  measured  bearings  in  the  interval  (t  ,  t  ) 
is  ^ 


mj+l+P 

■■j  - 

(C-1-3) 


in  which  the  B^^j's  are  measured  bearings  at  times  t^^  and  the  B*.  .'s 
are  predicted  bearings  for  times  t^.  ^ 

Two  schemes  have  been  considered  for  bearing  prediction.  In  the 
first,  the  predicted  bearings  are  obtained  from  the  quadratic 
expression 


(0-1-1) 


I 

I 

I 

I 

I 

I 

T 

I 

I 

I 


In  which  a j ,  bj,  and  Cj  are  obtained  from  a  least  square  fit  of  the 
m.  bearings  in  the  interval  (t  ,t.). 

°  .1 

Hie  random  part  of  the  deviation  r.  for  the  quadratic  extrapolation 

] 

:i 

(C-3-5)  n 


in  which  the  c^j's  are  as  given  in  equations  (A-1-12)  with  t^^  replaced 

I! 


OiiP^+CggT^jl+OjjT^  j2+2Ci2TjiP 

1/2 

■^^°13'^J2P+2°23'^Ji'^J2+P 
'  > 
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by  t  j ,  and  Tj^  and  Tjg 


(C-3-4) 


'The  alternate  scheme  Is  based  on  linear  prediction  from 


'tj  - 


(C-1-2) 


In  which  a'j  and  b'j  are  obtained  from  a  linear  least-square  fit  of 
the  mj  bearings.  The  Interval  considered  must,  of  course,  be  limited 
such  that  equation  (C-1-2)  is  a  sufficiently  accurate  description. 

The  standard  deviation  for  the  linear  extrapolation  Is 


(C-4-3) 


In  which 
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The  sums  In  the  above  are  from  to  tj.  The  linear  least-square 

fit  is  obtained  by 


a 

I 

:i 

] 

T 

1 

I 

1 

I 

I 


The  probability  of  a  target  zig  (P  )  is  found  by 


(C-1-4) 

3.  Theoretical  Results 

The  theoretical  results  of  Section  6  of  Appendix  C  Indicate  that,  at 
least  for  long  range  targets,  the  linear  prediction  scheme  Is  superior 
to  the  quadratic  scheme.  Further  study  at  the  shorter  ranges  Is 
needed  especially  where  the  time  for  executing  the  zlg  is  long.  The 
quadratic  scheme  might  be  better  In  such  cases.  It  Is  quite  possible, 
however,  that  less  sophisticated  schemes  such  as  visual  Inspection  of 
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a  bearing -time  curve  v/ould  be  completely  adequate  for  these  short 
ranges . 

Pigs.  C-6-1  through  C-6-5  (see  Section  6  of  Appendix  c)  and  the  ac¬ 
companying  tables  give  comparisons  of  and  Oj  for  various  geome¬ 
tries.  Approximate  expressions  for  andOj  v/ere  used  to  obtain 
these  figures.  The  theoretical  expressions  used  to  determine  are 
based  on  the  following  assiimptlons  : 

1.  Range  does  not  change  appreciably  during  the  time  Intervals 
considered. 

2.  Bearing-rate  for  each  leg  of  the  target  track  is  approxi¬ 
mately  constant. 

3.  Target  zlgs  occur  instantaneously  and  without  changes  In 
speed. 

Assumptions  1  and  2  above  are  not  unrealistic  for  long  ranges,  l.e., 
ten  nautical  miles  or  more,  itie  effects  of  assumption  3  can  be  par¬ 
tially  reduced  by  applying  the  following  Interpretations  for  "zlg 
time"  ! 

1.  The  "time  of  zlg"  (as  used  In  the  diagrams  and  tables)  should 
be  considered  to  lie  approximately  mldviay  between  the  begin¬ 
ning  of  the  zlg  and  the  end  of  the  zlg. 

2.  Speed  reductions  resulting  from  course  zlgs  advance  or  re¬ 
tard  the  "time  of  zlg"  depending  on  v^rhether  the  speed  reduc¬ 
tions  Increase  or  decrease  the  change  In  bearing-rate  caused 
by  the  change  In  course.  Thus,  apart  from  the  correction  In 
1  above,  "zlg  time"  should  be  Interpreted  as  the  time  at 
which  the  bearing -rate  changes  abruptly. 

The  values  of  and  consequently  the  zlg  probabilities  are  based  on 
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an  assumed  rav/  bearing  standard  deviation  of  0.2  degrees  and  a  bear¬ 
ing  sampling  Interval  of  2  seconds.  Larger  standard  deviations  v;lll 
degrade  the  results  somewhat  whereas  a  constant  bearing  off-set  will 
have  no  effect. 

The  geometries  considered  (Pigs.  C-6-1  through  C-6-5)  represent  zlgs 
which  are  particularly  difficult  to  detect.  In  all  cases,  a  15  de¬ 
gree  countercloc!<wlse  rotation  of  the  target  track  v/ould  make  the 
zlg  Impossible  to  detect.  In  such  cases,  hov/ever,  the  zlg  would  not 
affect  the  bearings -only  range  solutions.  Zlgs  v;hlch  cause  larger 
changes  In  bearing-rates  v;lll,  of  course,  result  In  larger  zlg  proba¬ 
bilities  at  a  given  time. 

Pigs.  C-6-1  and  C-6-2  and  the  corresponding  tables  Illustrate  the  su¬ 
periority  of  the  linear  zlg  detection  scheme  over  the  quadratic  scheme 
at  a  relatively  long  range.  Since  the  assumptions  used  for  determin¬ 
ing  r^.  are  not  valid  at  short  ranges,  a  similar  comparison  cannot  be 
made  for  the  short  range  cases;  computer  simulations  will  be  used  for 
this  latter  situation.  The  zlgs  for  these  tv;o  figures  v/ere  assumed 
to  occur  at  the  beginning  of  the  scan  Intervals. 

Pigs.  C-6-3  through  C-6-5  and  their  respective  tables  Illustrate  the 
growth  of  the  zlg  probability  with  subsequent  scans.  Thus,  If  the 
zlg  Is  not  detected  (zlg  probability  too  small)  In  the  Interval  In 
which  It  occurred.  It  will  be  easier  to  detect  In  the  subsequent  In¬ 
tervals.  This  will  not.  In  general,  be  true  for  the  quadratic  scheme. 

Analysis  Is  needed  to  determine  the  effects  of  undetected  target  zlgs 
on  the  range  estimate. 
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IV 

CONCLUSIONS  AND  AREAS  REQUIRING  FURTHER  INVESTIGATION 


1.  Conclusions 

/ 

A.  The  ASP  relative -motion  analysis  is  superior  to  manual  plot¬ 
ting  analysis  in  the  following  respects: 

1.  Accurate  relative -motion  parameters  are  obtained  in  less 
time  for  ASP  than  that  required  for  manual  plotting. 

2.  Ilie  ASP  system  can  analyze  targets  at  longer  ranges  than 
are  possible  using  manual  plots. 

B.  Itie  ASP  relative -mot  ion  analysis  is  superior  to  the  present 
Mark  113  system  analysis  in  the  following  respects: 

1.  More  relative -motion  parameters  are  obtained  as  outputs 
from  the  ASP  analysis  than  from  the  present  Mark  113  sys¬ 
tem  analysis. 

2.  The  ASP  system  has  the  capacity  to  produce  "if -therefore" 
solutions  based  upon  relative -motion  parameters  and  an 
hypothesized  input  of  target  range,  course  or  speed.  The 
Mark  113  system  can  obtain  this  type  of  solution  only  for 
an  input  of  target  speed.  Also,  the  possibility  of  a  dual 
solution  is  recognized  by  the  ASP  system,  but  Ignored  by 
the  present  Mark  113  system. 

C.  Automatic  target  zig  detection  is  possible  within  the  ASP  fire 
control  computer. 
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D.  The  linear  zlg  detection  technique  Is  superior  to  the  quadra¬ 
tic  technique  at  the  longer  ranges. 

E,  The  significance  of  curvature  criterion  for  establishing  the 
widths  of  bearing  pre-smoothlng  Intervals  Is  superior  to  the 
bearing-rate  criterion  presently  employed  In  the  Mark  113 
system. 


F.  As  presently  formulated,  the  ASP  bearlngs-only  complete  solu¬ 
tion  Is  Inferior  to  that  of  the  Mark  113  system. 


2.  Areas  Requiring  Further  investigation 

A.  A  criterion  Is  needed  whereby  the  Interval  over  which  the 
bearlng-tlme  expansion  Is  applied  Is  such  that  the  expansion 
la  Indeed  an  adequate  approximation  to  the  true  bearlng-tlme 
function.  In  the  event  that  the  total  Interval  can  only  be 
adequately  approximated  by  two  or  more  separate  expansions. 

It  would  be  desirable  to  combine  the  Information  obtained 
from  the  separate  Intervals  In  such  a  manner  that  no  useful 
Information  Is  discarded.  The  need  for  such  a  criterion  Is, 
as  previously  Indicated,  more  Important  at  the  shorter  ranges. 

B.  The  bearlngs-only  solution  as  presented  here  represents  a 
"first  look"  at  the  system.  The  asstunptlon  that  own-ship  zlg 
occurs  Instantaneously  must  be  replaced  by  a  realistic 
maneuver. 

C.  The  combination  of  the  first  and  second  leg  Information  has 
not  been  fully  exploited.  The  scheme  presented  here  treats 
the  six  parameters  (a's,  b*s,  and  c's)  as  Independent  whereas 
only  four  can  truly  be  Independent  (x,  y,  x,  and  y  of  the 
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target) .  It  would  be  desirable  to  find  relations  describing 
dependence  If  such  can  be  obtained  without  Introducing  the 
complexity  of  the  previous  SUBIC  schemes. 

Variance  expressions  for  many  of  the  parameters  have  yet  to 
be  derived. 

A  method  of  Introducing  additional  measurements  (e.g.,  range 
from  PUFFS  or  single  ping)  Into  the  bearings -only  complete 
solution  is  necessary  to  obtain  the  best  possible  solution 
from  all  of  the  data  which  might  be  available. 

Analysis  of  the  zlg  detection  schemes  must  be  continued  and 
extended.  Further  study  at  shorter  ranges  Is  needed,  es¬ 
pecially  where  the  zlg  Is  extended  In  time.  Analysis  Is  also 
needed  to  determine  the  effects  of  undetected  target  zlgs  on 
the  range  estimate. 

The  ASP  approach  to  fire  control  has  been  designed  to  permit 
more  effective  communication  between  the  man  and  his  computer. 
The  next  step  Is  to  Investigate  how  this  communication  can  be 
accomplished  and  to  evaluate  the  increase  of  tactical  effect¬ 
iveness  resulting  from  this  communication.  In  order  to  accom¬ 
plish  this,  it  is  planned  to  connect  a  research  console  to  a 
computer  to  determine  those  parameters  which  should  be  under 
operator  control  and  those  parameters  which  should  be  under 
computer  control. 
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APPENDIX  A 

DERIVATION  OF  TARGET  LOCALIZATION  EQUATIONS 


1 .  Least-Square  Approximation  to  the  Bearing-Time  Relation 

The  bearlng-tlme  relation. 


B 


B^+tan 

o 


(UAo^  % 

-'■-(u/RoHt:-to)cos 


(A-1-1) 


Is  obtained  from  Pig.  A-1-1  (which  for  convenience  Is  reproduced  here 
from  Pig.  1-1-2).  An  approximation  to  this  relation  Is  obtained  by 
taking  a  Taylor  expansion  of 


tan(B-B^) 

and 

(UAo)(t-to)s^ri 

-‘■-(U/RoHt-t^;cos 

to  get 

B  =  a+b(t-t^)+c(t-t^)^  (A-1-2) 

In  which 

^  (A-1-3) 

^  =  (U/RQ)sln  a^,  (A-1-4) 
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w  and 

c  =  (U/^?Q)^sln  cos  a^.  (A-1-5) 

When  the  quadratic  expansion  (equation  (A-1-2))  Is  sufficiently  accu¬ 
rate  and  when  the  bearings  have  random  variations,  a  least -square 
processing  of  the  bearing  equation  will  yield  unbiased  estimates  of 
the  coefficients  a,  b,  and  c. 

Assuming  that  the  error  In  time  Is  Insignificant,  the  sum  of  the 
squares  of  the  residuals  Is  expressed  by 


G 


a+b(tj^-t^)  +  c(t^-t^)‘^-B^ 


T  2 


(A-1-6) 


which  Is  minimized  with  respect  to  the  coefficients  to  give 


an+b2(tj^-tQ)+c2(t^-tQ)^  =*  SB^, 

a2(ti-to)+b2(ti-to)^+c2(ti-to)^  =  2B^(t^-t^),  (A-1-7) 

a2(tj^-t^)2+b2(t^-t^)3+c2(t^-t^)^  =  ^B^^^l-^o)^* 


Let  the  coefficients  of  a,  b,  and  c  be  designated  by  the  matrix  A, 
then 


a 

2Bi 

1«1 

b 

= 

2Bi(ti-to) 

c 

(A-1-8) 
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which  can  be  written  ae 


In  which 


(A-1-9) 


(A-1-10) 


V  = 


(A-1-11) 


Then 


C  =»  A 


(A-1-12) 


^  °21^1'^°22^2'^°23^3» 


Explicitly,  these  equations  are 


Det  |a  I  =  nI(t^-t^)2r(tj-t„)'‘+2t(t^-t„)£(tj^-t„)^E(tj^-t^)3+ 

-  ^-n 
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a  Det  I  A  (  = 

_  2 

b  Dst  |a|  -  nIBj^(t^-t^)r(t^-t^)V£B^(tj-t^)2l(ti-t^)t(tj-t^)2+ 

r  12 

-n£B^(  tj^-t^)2£  ( tj-t^)3.£B^£(  t^-t^)£  ( t^-t^)  \ 

and 

c  Dot  |a|  =  nlB^(t^-t^)^Z(t^-t^)^+ZBj^Z(t^-t^)Z(tj^-t^)3+ 

r  1 2 

+  ZBi(ti-to)Z(t^-to)l(ti-to)^-EBi[z(t^-t^)2|  + 

r  1 2 

-n£Bj^ ( t -t^ )  £(  )  3-  £Bi ( t  j -t^ )  2  [£  ( )  . 

2.  Beai*fng-Rate  and  Change  In  Bearing-Rote 

.  •:a  vin.^-rate  ,  the  time -rate -of -clianpe  of  bcarlnj^,  is  found  by  dlffer- 
^r■,  eo'-.,?  t  i.(,n  (A-.t-i')  v.'lth  respect  to  time.  Thus, 
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similarly,  the  change  In  bearing -rate  Is  given  by 


B  =  2c. 


(A-2-2) 


3.  Determination  ofPT.  9  and  U/R^ 

Prom  equations  (A-1-4)  and  (A-1-5), 


a  tan~^(b^/c) 


(A-3-1) 


and  from  Fig.  A -1-2, 


a(t)  =  a^+B(t)-a. 


(A-3-2) 


This  angle  is  referred  to  as  the  relative  angle -on-the -bow. 
The  relative  course  (9),  see  Pig.  A-1-1,  is  seen  to  be 


From  equation  (a-1-4). 


0  =  TT  -  -f  a. 


U/Tl^  =  b/sln 


and  from  equations ( A-3-1)  and  (a-3-4). 


(A-3-3) 


(A-3-4) 


2 

c  +b 


(A-3-5) 
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4.  Development  of  Solution  Using  Speed  Input 

data,  own-ship  motion  parameters  and  a  given  target  speed 
(Vg'  /  used  to  determine  a  complete  solution. 

Prom  Pig.  A-4-1  (which  Is  Pig.  1-1-3  reproduced  here  for  convenience) 
target  course  Is  given  by 


-  T-r+a-a 

o 


(A-4-1) 


In  which 


sin 


-1 


V 


(1) 


Je 


m 


sin 


a+“o) 


(A-4-2) 


The  foregoing  equation  Indicates  the  possibility  of  two  solutions 
(this  Is  a  property  of  all  target  localization  solutions  based  on 
bearing  data).  Pig.  A-4-2  Illustrates  the  fact  that  when  bearing 
lines  converge,  two  possible  target  tracks  exist  for  a  given  target 
speed.  By  examining  the  angle  (jZf^ -a+a^)  the  various  solution  pos¬ 
sibilities  may  be  determined. 


Pig.  A-4-3  (a)  and  (b)  Is  a  vector  velocity  diagram  In  which  the 

angle  (;zr^  '-a+a^)  Is  In  the  first  and  fourth  quadrants,  respectively. 
It  can  be  seen  that 


min  '' 


sln(;zr^^)-a+a  ) 


(A-4-3) 


When  exceeds  but  is  less  than  v^^),  two  values  of  r  exist. 
One  Is  given  by  the  principal  value  of  equation  (A-4-2)  and  the  other 
by  Its  complement.  TSius, 
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%-«#•  -a  +  aQ^2Tr 

Fig.  A-4-3  First  and  Fourth  Quadrant  Solution 
Possibilities 
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>1  =  r 


(A-4-4) 


^2  "  TT-y, 


(A-4-5) 


When  only  exists  and  there  Is  a  unique  solution. 

Fig.  A-4-4  (a)  and  (b)  depict  the  situations  In  which  the  angle 
(^(^)_a+a^)  lies  within  the  second  or  third  quadrants,  respectively, 
The  minimum  target  speed  Is  given  by 


min  • 


(A-4-6) 


For  any  larger  value  of  a  unique  solution  exists  and  (prin¬ 

cipal  value  of  equation  (a-4-2))  alone  exists. 

Fig.  A-4-5  sxommarlzes  the  various  solution  possibilities. 

When  two  values  of  7  exist,  the  entire  solution  must  be  performed  In 
duplicate,  llius,  equation  (A-4-1)  becomes 


^(3) 


(J  »  1  or  1,  2).  (A-4-7) 


From  the  law  of  cosines 


Uj  =  (A-4-f 


and,  from  equation  (A-3-5) 


P  li  1/2 

Rqj  =•  bUj/(c  +b  ) 


(A-4-9) 
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Fig.  A-4-5  Speed  Input  Solution  Possibilities 
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For  range  at  any  time  (t).  It  can  be  seen  from  Pig.  A-1-1  that 

U  sin  cto(t-t^)  =  R(t)3ln[B(t)-a] . 

Thus,  assuming  (B-a)  to  be  small 


R(t)j 


R  .b 
oj 

b+c(t-t^) 


(A-4-10) 


It  la  possible  that  one  of  the  two  solutions  will  be  physically  ab¬ 
surd,  although  mathematically  possible. 

5.  Development  of  Solution  Using  Range  Input 

At  any  time  (t^),  consider  a  coordinate  system  with  the  x  axis  perpen¬ 
dicular  to  the  bearing  line  and  the  y  axis  along  the  bearing  line. 

The  components  of  target  speed  at  this  particular  time  can  then  be 
written  In  terms  of  knovm  parameters  and  a  range  input  at  the  given 
time,  Rg(tg).  The  cross-llne-of -sight  component  of  target  speed  at 
tg  Is  then  given  by 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

J 

* 

.1 

.1 


In  which 


B(tg)  =  a-f-b(tg-t^)+c(t^-t^)2  (A-5-2).  ■] 

and 

n 

B(tg)  =  b+2c(tg-t^)  .  (A-5-3) 

The  along-llne -of -sight  component  of  target  speed  Is  given  by 

li 
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=  V^^^cos[j^(^^-B(tg)]+Rg(tg)  (A-5-^0 


In  which  the  quantity  RgCt^)  can  be  developed  as  follows. 
Prom  equation  (A- ■^-10) 


Rob 


ITius , 


R^bc 

^e^^e^  -  7r~T7 — TTiiT  ' 


or 


«e(^e) 


[b+c(tg-to)l‘ 


(A-5-5) 


Equations  (A-5-1)  and  (A-5-^)  can  then  be  v;rltten  as 


-  v‘^>3ln[^(l)-B(t^)]+R^(t^)[tH-2o(t^-t^)]  (a-5-6) 


and 


Vy'^'(t3)  =  v(^’oos[)i((l)-B(t3)]-  .  (A-5-7) 


e  O' 


Total  target  speed  Is  found  from  Its  components. 


V 


(2)  _ 


X 


1/2 


(A-5-8) 
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From  Pig.  A-^-1,  It  can  be  seen  that  target  course  is  given  by 


=  B(t  )+tan'‘ 


iK 


(te) 


y  ^  e ' 


(A-5-9) 


Prom  equation  (A-^-lO) 


and  for  any  time  (t)  other  than  (t^) 


R(t) 


■nius,  range  at  any  time  is  given  by 


R(t)  = 


(A-5-10) 


6.  Development  of  Solution  Using  Course  Input 

data,  own-ship  motion  parameters  and  a  given  target  course 
i0Q^  )  are  used  to  determine  a  complete  solution. 


Prom  equation  (A-4-1) 


y  =  TT+a-a  -0  . 


(A-6-1) 


Itien,  rearranging  equation  (a-4-2) 
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V 


(2)  V^^^3ln(^(^)-a+ap) 

sin  y 


(A-6-2) 


Target  range  can  then  be  obtained  from  equation  (A-4-10)  In  which 


(a-6-3) 


7.  Development  of  Bearings-Only  Solution 

Own-ship  travels  and  collects  bearing  data  on  two  distinct  "legs", 
each  consisting  of  uniform  motion,  but  differing  In  speed  and/or 
course.  The  simplifying  assumption  Is  made  that  own-ship  zlg  (change 
of  course  and/or  speed)  occurs  Instantaneously,  The  subscripts  1  and 
2  Indicate  first  and  second  leg  parameters,  respectively.  The  x  and 
y  subscripts  Indicate  the  cross-llne-of-slght  and  along-llne-of -sight 
speeds,  respectively.  iTie  Initial  time  Is  designated  by  t^,  the  time 
of  own-ship  zlg  by  t',  and  current  time  on  the  second  leg  by  t. 

At  the  time  of  own-ship  zlg,  the  relation  existing  between  range, 
relative  cross-llne-of-slght  speeds  and  bearing-rate  Is 

R(t')B(t')  =  U^(t')  =  V^^^)(t')-V^^^^t').  (A-7-1) 


For  a  uniform  target  track. 


R(t')AB(t')  =  -  AV^^^^(t'.). 


(A-7-2) 
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Thus, 


R(t') 


B2{f  )-Bi(f  ) 


(A-7-3) 


In  which 


Vax'^’C:')  -  V2'^'3ln[)i(2(^)-B(t')],  (A-7-5) 


and  the  bearing  at  time  t'  la  the  time-weighted  mean  of  the  two  bear¬ 
ings  (one  from  each  of  own-ahlp'a  "legs")  obtained  at  time  t'.  Thus, 


B(t')  =  |-^[ai+b3^(t'-tQ)+c^(t'-t^)^]  + 


t-t' 


(A-7-6) 


The  two  bearing-rates  at  time  t'  are  given  by 


B^(t')  =  b^+2o^(t ' -t^) 


B2(t')  -  b2. 


(A-7-7) 


(A-7-8) 


Range  at  any  time  on  the  second  leg  (t)  Is  then  calculated  from 


R(t) 


R(t')b2 


(A-7-9) 
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. : 


By  assumption, targe t  course  and  speed  are  constant.  Therefore,  In¬ 
formation  obtained  on  each  leg  of  own-ship  track  can  be  used  to  cal¬ 
culate  target  course  and  speed  by  the  method  used  In  the  range  Input 
solution.  Thus,  for  leg  k  (k  =  1,  2), 


)  =  Vj^^^^sln 


^j^(^^-B(t')  +R(t')^(t')  (A-7-10) 


and 


Vky^^^(t')  =  Vj^^^^cos  4^^)-B(t')]+Ri,(t' )  (A-7-11) 


In  which  B(t')  Is  given  by  equation  (A-7-6), 


R^(t') 


R(t')c^ 

fel+Ci(t' -to)  " 


and 


R(t')c2 

«2(^')  -  - 


(A-7-12) 


(A-7-13) 


Target  speed  Is  then  given  by 


and  target  course  Is  given  by 


=  B(  t '  )  +  tan"^ 


(A-7-15) 
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The  results  of  equations  (a-7-1^)  and  (a-7-15)  are  then  weighed 
to  obtain  the  best  values  of  target  speed  and  course.  Thus, 


-  WjVj(2)  +  (A-7-16) 

and 

0^^'^  =  ’^3^1^^^  +  (A-7-17) 


In  which  the  W's  are  the  appropriate  weighing  factors. 


1 

:i 

] 

I 


1 

! 

..  & 


! 
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APPENDIX  B 

DERIVATION  OF  STATISTICAL  PROPERTIES 


1 .  Coefficient  Variances  and  Covariances 
By  definition. 


C  =  AdJ  A/Det  )  A  | 


(B-1-1) 


in  which  A  la  the  coefficient  matrix  of  equation  (A-1-8)  and  let  c. - 
th  ~  iJ 

be  the  IJ  element  of  C.  Since  A  =  A,  C  =  C,  i.e.,  both  matrices  are 

symmetric,  or  c^j  =  then. 


variances 


°a  “  °11°B 


°b  “  °22°B 


°c  “  °33°B 


(B-1-2) 


covariances 


‘'ab°a°b  "  °12°B 


■'bc^b^c  “  °23°B 


T  O  O  =  p  O  ^ 

ac  a  c  ^13  B 


(B-l-3) 


2 

in  which  Og  is  the  variance  in  the  bearing  data  and  the  are  the 

correlation  coefficients  which,  from  equations  (B-1-2)  and  (B-l-3)  are 


63 


CONFIDENTIAL 


CONFIDENTIAL 


■^ab  “ 


(C22C33) 


r,  and  T 


°  7o,^c33)Vi;- 


mie  elements  of  A  from  whleh  the  c^j-e  are  oaloulated  consist  of  rel¬ 
atively  oumbersosie  sumnatlone.  Simplified  approximations  are  presented 

here  to  facilitate  calculations.  Let  the  bearings  be  obtained  at  equal 
time  Increments  (T),  then 


(B-1-5) 


in  which  t^  is  the  time  at  which  the  bearing  is  obtained. 


Also, 


Similarly, 


“  «  n 

-  rTl»Tfl 
1=“0  l*«o  Iwo 


£t^  -  T  n/2  (n+1) . 


1^1^  -  n/6  (n+l)(2n+l), 

It^3  .  .J.3  (n+l)2, 

It^^  -  nVlO  (n+l)2(2n+l) 


(B-1-6) 


For  n  large , 


Lt^  =  T  n/2. 


nVs, 


It^^  =  n^/k,  and  It^^  *  n^/5. 


(B-1-7) 
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Relating  equations  (B-1-7)  and  (a-1-7),  the  elements  (a. .)  of  A  are 

J-  J 


^12  ^21 

^13  ^  ^31  “  ^22  “  ^  (B-1-8) 

O  ll 

^2^  ~  ^32  ^  T  n  >  3.nd 

^33 " 


The  determinant  of  A  is 


Det  I  A  1  =  +  a^^a^^a^^  +  a^^a^^a^^  + 

-  S-^1^22^13  ■  ^32^23*^  ■  ^21^12^33 

and  using  (B-1-8), 

Det  |a  I  =1  nV{2l60).  (B-1-9) 

The  elements  of  C  are. 


11  ^  ^^22^33 

■*  a^2e2^)/Det 

1  a|. 

22  “  ^^11^33 

~  ^31^13)/^^^ 

1  a|. 

(B-1-10) 

33  “  ^^11^22 

1  a|. 
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I 

1 

C21  =•  0^.2  “  I^I*  -| 

®32  ”  °23  “  (®'31®^12”^11®'32^^^^^  (B-1-10)  -r 

and  J 

®13  "  °31  “  (^12^23“^22^13^/^^  |a|.  -i 


Using  equations  (B-1-8)  and  (B-l-9)» 

;i 

°11  “ 

] 

°22  “  192/T^  n3> 

0^3  -  180/T^  n^. 

(B-1-11) 

1 

-T 

C21  “  c^2  “  -36/T  n  , 

] 

^23  *  ^32  "  “180/T^  n  f 

] 

and  “  °3i  “  30/T^  n^. 

1 

Combining  equations  (B-1-2),  (b-1-3),  and  (B-1-11), 

1 

Oa^  -  (9/n)Og2, 

:i 

0/  -  (192/T^n3)Og2^ 

(B-1-12) 

'1 

0/  -  {l80/TV)Og2^  _ 

j 

i 

and 

1 

’abVb  " (-36/1^2) "b®- 

;! 

’boS'c  =(-l8o/T3n'*)o^2^ 

'ao^a'c  =(30/T®n3)<.g2_ 

(B-l-13) 
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Also,  from  equations  (B-1-4)  and  (B-l-ll), 


-VV2  =  ~  0.866, 


-^^T5/h  «  -  0.970, 

VT/s  -  0.745. 


(B-1-14) 


Pigs.  B-1-1,  B-1-2,  and  B-1-3  are  plots  of  the  standard  deviations  of 
the  least  square  coefficients  as  functions  of  time.  A  bearing  sampling 
interval  of  two  seconds  and  a  sonar  bearing  standard  deviation  of  0.2 

degrees  are  used. 


2.  Variances  In  Functions  of  the  Coefficients 

For  a  linear  combination  of  the  coefficients 


u  =  k^a  +  kgb  +  k^c, 

“  '^1  \  +  2^ab‘^l''2‘’a°b  + 


(B-2-1) 


For  an  arbitrary  function  of  the  coefficients 


u  =  f(a,  b,  c). 


=  (^)2o  2  /^\2  2  2  2  Su  5u 


I  o  *T  ^  ^Vl 

\r.  TTh  5F  °b®c  ^ 


ac  5a  5c 


(B-2-2) 
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to  a  first  order  approximation. 

Variance  expressions  for  some  of  the  functions  of  Interest  are  devel¬ 
oped  In  the  following  pages.  Analysis  Is  being  continued  on  several 
that  have  not  yet  been  evaluated. 

3.  Variance  In  Bearing  Rote  and  Chcmge  In  Bearing  Rote 


Prom  equation  (A-2-1), 


B  ■  b  +  2ct. 


Iberefore,  taking  =■  -1  (equation  B-1-14), 


SF  "b°c 


„  dB  „ 

°B  SF  ®b  "  5F  0  » 


which  results  In 


“b  -  %  -  • 


(B-3-1) 


Prom  equations  (B-1-12), 


°b  "  ^'^c' 


or,  from  the  definition  of  T  and  n. 


%  =  ■^°c- 


71 


CONFIDENTIAL 


CONFIDENTIAL 


Therefore, 


°B  “  ‘^b-  (B-3-2) 

Hius  Pig.  B-1-2  can  be  Interpreted  as  a  plot  of  versus  tracking 
time.  Ihe  time -rate -of -change  of  bearing -rate  Is  given  by 

B  -  2c.  (A-2-2) 


QSius,  It  Is  obvious  that 


o 


*  • 

B 


(B-3-3) 


Fig.  B-3-1  Is  a  plot  of  og  versus  tracking  time. 


4.  Varlancot  ln<i^Q,o(,  and  9 

Prom  equation  (A-3-1) 


=  tan"^(bVc) 


and 


o  2b/ c _  2bc 

^  [i+(bVc  )^]°  ■ 

_  -b^/e^  .b2 

^  [Cb^/c)^] 

Using  equations  (A-l-^t)  and  (a-1-5) 
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=  2(u/k^)“-^cos  and 


“  -  (uAo)"^- 


Thus,  using  =  -1, 


2(U/R^)“^cos  a  o  +(u/il^)“2o 


O'  c 


(B-4-1) 


Consider  the  ratio  of  the  first  term  in  (B-4-1)  to  the  second  term. 


2(UA^)"-^cos  a  o, 


V  I  uua  Ok 

®  ^  _  O  /  TT  As  ~  D 


(U/fo)  “c 


=  2(U/^  )co3  a  —  . 
'  O  o  o 

c 


Using  (B-l“12),  the  ratio  becomes 


2(U/^q)cos  a^(Tn)  =  2{\J/R^)t  cos 


which  must  be  small  for  the  expansion  expressed  by  (A-1-2)  to  be 
valid.  Thus,  to  a  first  order  approximation. 


-a  "  (UAo)  ^ 


(B-4-2) 


It  can  be  shown  that 


o  «a  o  . 
cz  a 

o 


(B-4-3) 


Prom  equations  (a-3-3)  and  (A-3-1) 


e  «=>  ir-tan“^(b^/c)+a. 
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Also, 


and 


he 

3a 


1, 


he 

3F 


de 

3c 


-2bc 

• 


using  -  -1,  =  -1,  and  -  +1, 


,  2bc  b*^ 

^  ^2..  4  b  ^2,^  c 

c  +b  c  +b 


or 


‘’a  +  TTT  (2cOb+bo  ) 

c  +b 


(B-4-4) 


Consider  the  ratio  of  the  first  term  In  the  above  parenthesis  to  the 
second  term,  using  equations  (A-1-4),  (a-1-5),  and  (B-1-12), 


2 

2co^  2(U/^^)  sin  cos  to^ 

"  tF/h„  sin 

C  '  o  o  c 


2(U/^„)cos  a  t 
'  '  o'  o 


which  must  be  small  for  the  bearing  expansion  to  be  valid.  Itierefore, 
to  a  first  order  approximation. 


o 


e 


■  2^.  4 

c  +b 


(B-4-5) 
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TO  Versus  Time 
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Then,  from  equation  (A-3-5), 


_  /R  \2 

O.  a  O  +  o 

6  a  U  c 


(B-4-6) 


Expressing  in  terms  of  and  using  the  arg\aments  used  above 


o  ^ 

®  1^0 


o  =  o  . 
c  a 


(B-4-7) 


Pig.  B-4-1  Is  a  plot  of  and  versus  tracking  time  for  various 
values  of  U/R^ . 


5.  Variance  In  (U/R^) 

Prom  equation  (A-3-5) 


c  +b 


Using  equations  (A-1-4,  5)  and  the  above. 


S(UAq) 


.  2  2 
sin  a  -cos  a 
_ o _ o 

sin  a 


^(U/Rq) 


cos  a 


^  =  (U/R^)-^cot  a 


Therefore,  taking  =  -1, 
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Pigs.  B-5-1  through  B-5-3  show  the  variation  of  Oy/Ro  tracking 

time  for  various  values  of  and  a^. 

For  =  Tr/2, 


j  2  2  \ 

sin  a^-cos 

sin  , 


-1 


cot  0.0 
o  c 


(B-5-1) 


(B-5-2) 


For  0^=0, 


a 


oo  , 


(B-5-3) 


In  this  case,  however,  bearing-rate  Is  zero  and  the  only  Information 
which  can  be  derived  from  the  data  Is, 


B 


constant  and  B  =  0. 


6.  Examples  for  Particular  Geometries  and  Comparisons  of  ASP  with  Manual  Plots 

The  figures  of  this  section  give  relative  motion  parameter  accuracies 
as  functions  of  time  for  the  geometries  Illustrated.  The  ASP  results 
are  based  upon  a  two-second  bearing  sampling  Interval  with  a  bearing 
standard  deviation  of  0.2  degrees,  and  the  manual-plot  results  are 
based  upon  a  one  minute  bearing  sampling  Interval  with  a  bearing  stan¬ 
dard  deviation  of  0.3  degrees.  The  significance  of  these  figures  Is 
discussed  In  Section  ^-A  of  Part  II. 
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Fig,  6-6-1  Case  1  Geometry 
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Fig.  B-6-3  Case  1  cr^  Versus  Time 
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CASE  2 
GEOMETRY 


Fig.  B-6-6  Case  2  Geometry 
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CASE  3 
GEOMETRY 


Fig.  B-6-1 1  Cate  3  Geometry 


TARGET 
10  KNOTS 
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CASE  3 


IN  DEGREES 


Fig.  B-6-.13  Cases  tr^Versus  Time 
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Fig.  B-6-16  Case  4  Geometry 
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CASE  4 

VS.  TRACKING  TIME  IN  MINU1 
Oq  =  I09* 


TRACKING  TIME  IN  MINUTES 


ig.  B-6-18  Case  4  (Tq Versus  Time 
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IN  PERCENT 


Fig.  B-6-20  Case  4 Versus  Time 


101 


CONFIDENTIAL 


CONFIDENTIAL 


I 


TARGET 
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CASE  5 
GEOMETRY 
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Fig.  B-6-22  Case  5  Time 
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10‘5 


IN  PERCENT 


TRACKING  TIME  IN  MINUTES 

Fig.  B-6-25  Case  5  Versus  Time  | 
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7.  Bearings-Only  Range  Variance 

From  equations  (A-7-3,  7,  8,  9), 


AV  )bp 

R(t)  - - - = - ^  . 

[b^+2c^(t' -tQ)-b2] [b2+C2(t-t ' ) ] 

Prom  this  equation. 


Sr{  t) 

-  ,  . 

R{t) 

R(t) 

6b2 

-  b,_,  + 

[b^+2c^(t'- 

^o)-^2]  [b2+C2(t-t')] 

aR(  t) 
dC2 

R(t)( 
02+02 ( 

:t-f)  » 

SR(t) 

-  R(c) 

_  R(t) 

8b^ 

[bi+2c^(f-t^)-b2) 

AB(t') 

dR(  t) 

-2R(t)(t '-to) 

2R(t)(t'-t^) 

dc^ 

^^3^+2c2^(  t  ‘  “t^)  -b2  j 

AB(t') 

Own-ship  speed  and  course  are  assumed  to  have  negligible  error  and  the 
own-ship  track  Is  assumed  to  be  nearly  perpendicular  to  the  B(t')  bear¬ 
ing  line  such  that 


VjU)sln(Aj+ej,) 


107 


CONFIDENTIAL 


CONFIDENTIAL 


Vj'  '^sin  Aj+Eg  cos  Aj J 


which,  to  a  first  order  approximation,  is 


Aj. 


Tlius,  errors  arising  from  resolving  ovm-shlp  velocity  are  negligible. 

Prom  the  above,  using  =  ^b2C2  ""  correlation  between  first 

and  second  leg  coefficients  and  =  to^. 


n  2  ^2 

R(t)^  AB^(t') 


in  which  Aj  *=  Tr/2  and  Cg  represents  the  error  in  B(t').  Thus, 


and  using 


A  B  =  B2  -  B^j^ , 


Bg  “  b2 } 


°R(t)  _  1  L  2 


2  .lA 


(B-7-2) 


Using  (B-1-12), 
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(B-7-3) 


in  which  and  Tg  are  the  times  on  the  first  and  second  legs,  respec 

tlvely,  and  T  is  the  bearing  sampling  Interval.  Let  the  bracketed  ex 

2 

presslon  be  denoted  by  z  ,  then 


(B-7-4) 


This  equation  is  similar  to  that  for  range  variance  in  Mode  2  (Refer¬ 
ence  l).  The  optimum  value  of  z,  however,  differs  In  the  two  systems. 

The  tracking  time  on  each  leg  appears  only  In  z.  Now  let 


Tg 

X  =.  =-  and  y 
■^t 


B. 


®l"®2 


In  which  Is  the  total  tracking  time.  Then, 


(1-x) 


1  .1 
I  3 


(y-i)‘ 


Minimizing  this  expression  with  respect  to  x  gives  the  optimum  time 
ratio  (Tg/T^)  and  the  minimum  z  (z-opt).  Figs.  (B-7-1)  and  (B-7-2) 
are  graphs  of  x-opt  and  z-opt  versus  y  for  both  this  system  and  Mode 
2  (Mbrk  113).  Z-opt  and  consequently  the  range  error  for  ASP  Is  about 
1 . 4  times  larger  than  the  error  for  Mode  2 . 
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ASYMPTOTE  AT  Z  •  4.00 


J'^ONFIDENTIAL 


111 


CONFIDENTIAL 


Fig.  B-7-2  z-opt  Versus  y 


CONFIDENTIAL 


Prom  equation  (B-7-2),  It  la  apparent  that  If  the  second  leg  bearing- 
rate  approaches  zero,  the  range  deviation  Increases  without  bound. 
Thus,  R{t)  may  be  of  no  value.  However,  the  range  at  time  of  own^ship 
zlg,  R{t'),  can  be  used  In  this  situation. 

From  equations  (A-7-3),  (A-7-7),  and  {A-7-8)  range  at  time  of  own-ship 
zlg  is  given  by 


B2(t' )-B3^(t' )  ^2“^l“^°l(^'-'‘=o^ 


Prom  the  above  equation. 


SR(t' ) 


5R(t' 


[b2-b^-2oj(t'-t„)]^  [AB(t')] 

[ b2-b^-2c^( t ' -t^) I ^ 


R(t' 


[AB(t')] 


5R(t')  _  ^(^'-^o)[Vix^^^(t')-V2^^^^t')]  _  2(f-t^)R{f) 
°1  [b2-b^-2c^(t'-t^)]^  [AB{t')] 

The  assumptions  which  led  to  (B-7-2)  are  again  used  to  get 
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2 

°  R(t 


f  R(t' )  1 

2 

» 

f-R(t' )  1 

2 

LAB(t')  ^1. 

+ 

LAB(t‘)  2. 

AB(t')  °°1_ 

T  2 


R(t') 

UB(t')J 


bi 


Thus, 


°%(t' ) 


R  Jt«) 


2  2 
o 


C  o 


Then,  using 


or 


»bj  -  (‘'-‘oJ'ci 

2  2 

q^R(t ' )  ^  °  bi'*'®  b2 
R^(t')  AB^(t') 


2  2 

1/2 

1  ^2J 

AB(t') 


Using  equations  (B-l-i2) , 


(B-7-5) 


(B-7-6) 
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APPENDIX  C 

DERIVATION  AND  ANALYSIS  OF  2IG  DETECTION  EQUATIONS 


1 .  Predicted  Versus  Measured  Bearings 

Let  the  Interval  t^  to  t^  be  such  that  a  quadratic  approximation  to 
the  bearlng-tlme  curve  Is  sufficiently  accurate,  and  let  the  bearing 
data  from  t^  to  tj  (t^  <  t^  <  t^)  be  least-square  fitted  to  obtain 
coefficients  Sj,  bj,  and  Cj.  The  bearing  at  time  t^  can  then  be  ob¬ 
tained  from 


(0-a-i) 


The  B  ^j's  thus  represent  predicted  bearings  based  on  bearing  history 
from  t^  to  ty  When  the  Interval  t^  to  t^  Is  further  restricted, 
bearings  can  be  predicted  from 


a 


+  b' 


(ti-to) 


(C-1-2) 


In  which  a'j  and  b'j  are  obtained  from  a  linear  least-square  fit. 
Let 


r 


J 


m.+p+l 

I 

l-eij+l 


-®lj) 


(C-1-3) 


In  which  the  B^^j's  are  the  bearing  data  as  received  from  sonar  at  t 
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Curves  for  a  Zigging 


ONFIDENTIAL 


nij  Is  the  number  of  bearings  prior  to  tj,  and  p  Is  the  number  of  pre¬ 
dicted  bearings.  It  can  be  seen  that  Tj  represents  a  quantity  propor¬ 
tional  to  the  area  between  the  expected  and  the  true  bearlng-tlme  curve 
from  tj  to  (see  Pig.  C-1-1).  This  area  Is  also  approximately 

proportional  to  the  difference  between  the  expected  and  the  true  bear¬ 
ing  rates. 

The  probability  of  a  target  zlg  Is  found  by 


(C-1-4) 


In  which  Oj  Is  the  standard  deviation  of  and  remains  to  be  deter¬ 
mined.  Zlg  detection  sensitivity  will  obviously  depend  upon  r  and  o.. 

J  J 

2.  Linear  Least-Square  Fit 

The  primed  coefficients  of  equation  (C-1-2)  are  obtained  from 

(C-2-1) 

(C-2-2) 


V' 


(C-2-3) 


(C-2-4) 
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(C-2-5) 


The  primes  are  used  here  to  distinguish  between  the  linear  and  quad¬ 
ratic  least-square  fits. 


Prom  the  above 


a*  =  +  °'l2^1^^1‘'‘^o^ ' 


b'  «■  ^  ®  *  22^®1  ^  ^l”^o^ 


(C-2-6) 


2  .  «  2  , 
-  °’ll°B 


2  2 
°b'  “  °'22°B  ' 


(C-2-7) 


■^a'b'°a'°b'  “  °'l2°B 


Prom  equations  (C-2-4,  5) , 


Det  A*  -  nI(tj^-t^)2-[E(t^-t^)]  , 


'll  -  ^{t^-t^r/Det  A'  , 
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and 


c'22  “  n/Det  A* 


(C-2-8) 


12 


21 


-Z(t^-t^)/Det  A' . 


1 

I 

1 

I 

} 

i 

1 


I 

1 

r 

r 

i: 

r 


3.  Variance  of  T  for  Quadratic  Extrapoiotion 

Prom  equations  (C-1-1,  3) 


(C-3-1) 


or 


ajP+bjl(t,.to)^CjE(t,-to) 


(C-3-2) 


In  which  the  sums  are  from  t^^  »  tj  to  Since  the  bearings 
are  bearing  data  obtained  after  time  tj  and  the  coefficients  Sj,  bj; 
and  c.  are  derived  from  bearing  data  obtained  prior  to  t,,  the  B. 


.j  - - - - - - - j, 

bearings  are  stochastically  Independent  of  the  coefficients.  Also, 

2 

Tj  Is  a  linear  combination  of  variates.  Thus,  the  variance  of  ^j(Oj  ) 
Is 


2  2 
/  P  '’a  + 


Z(ti-to) 


«  2 
%  + 


■  2 

‘’c  +  N 


{  +2T^^pL(t^-t^)o^o^+2T^^pZ(t^-t^)2o^o^  +  S.  (C-3-3) 


^  +2T^,^E(t^-t^)£(t^-t^)2o^o^+po 


B 
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Now  let 


I  (ti-t„)  - 


(C-3-^) 


ti=t^+T 


ti-tj-HT 


'^J2' 


then,  using  equations  (B-1-2,  3)  with  t^^  replaced  by  ty  equation 
(C-3-3)  becomes 


2  2 


CiiP  +C22'*’jl  ■^°33'^J2^'^^°12'^JiP 


■^^®13'^J2P'^^°23'^Jl'^j2'*'P 


.  (C-3-5) 


Let  the  number  of  bearings  from  t^  to  tj  be  m^,  then  the  of 

equation  (C-3-5)  can  be  replaced  by  equations  (B-l-ll)  In  which  n-m, 
For  t  -0  the  expressions  (C-3-^)  can  be  written  as 


nij+l+p 

Tjl  -  T  Y.  1 


mj+l+p 
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r 

X 

T 

* 

f 

t 


I 

L 


i; 

r 


r 

r 

I 


Now  let 


l=»K+nij+l , 


then  for  nij  large. 


and 


I  (K+mj)  =  Tfmj 
K=*0 


p+EK 

K 


P 

(K+m. 

K=«0 


)^=T^ 


m .  p+2nij 


ek+lk^ 

K  K  -I 


Ihe  stima  over  K  can  be  reduced  by 

.  P 

y  K  -  p(p+i)/2 

K-0 

and 


P 

=  p(p+l)(2p+l)/6. 

K-0 

Using  the  above  simplifications,  equation  (C-3-5)  becomes 


2  2 

‘’j  “‘’b 


2 


>  .  (C-3-6) 


i  ^^(p+l)^(2p+l)+  ^  (p-l)^(2p+l)^ 

"’J 

This  equation  was  used  to  obtain  the  curves  of  Figure  C-3-1. 
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STANDARD  DEVIATION  P  IN  DEGREES 


Fig.  C-3-1  Standard  Deviation  of  T  Versus  Time 

for  Quadratic  Extrapolation  CONFIDENTIAL 
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STANDARD  DEVIATION  OF  T  VS.  TIME  FOR 
_ LINEAR  EXTRAPOLATION 


Fig.  C-4-1  Standard  Deviation  of  T  Versus  Time 
for  Linear  Extrapolation 


123 


CONFIDENTIAL 


CONFIDENTIAL 


4.  Variance  of  Fj  for  Linear  Extrapoiotion 

Prom  equations  (C-1-2,  3) 


•  1 
I 


{°'iiP^'*'°'22'^j/+2c'^2TjiP+p}  (C-4-3) 


In  which  Tj^  la  defined  by  equation  (c-3-4) . 

Ttie  approximation  to  (C-4-3)  analogous  to  (C-3-6)  Is 


°/*  ”b^  {p*  ST  ^  (P+l)+  ^  (P+I)^  i  .  (0-11-4) 

J  nij  > 

Curves  of  Fig.  (C-4-1)  were  obtained  from  equation  (C-4-4). 


5.  Approxinnote  Expressions  forF^ 

llie  analysis  of  this  section  Is  simplified  by  assuming  target  zlgs  to 
occur  Instantaneously  and  without  changes  In  speed  for  course  zlgs. 
The  effects  of  these  assumptions  must  be  considered  when  Interpreting 
the  results  of  the  analysis.  Thus,  a  finite  zlg  time  will  retard 
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detection  while  a  speed  reduction  accompanying  a  course  zig  will 
either  retard  or  advance  detection  depending  upon  its  effect  on  the 
change  in  bearing-rate.  The  effects  can  easily  be  deduced  for  the 
particular  cases  examined. 

Consider  first  the  linear  extrapolation  scheme  expressed  by  equation 
.  This  scheme  will  obviously  be  suitable  for  long  range  tar¬ 
gets,  but  can  also  be  used  at  the  shorter  ranges  by  restricting  the 
time  intervals.  A  target  zlg  is  Illustrated  by  the  bearing-time  curve 
of  Pig.  C-5-1.  The  curve  consists  of  two  straight  lines  (lines  1  and 
2)  meeting  at  (the  time  of  target  zig)  .  The  case  where  (t^.  <  tj) 
will  be  considered  later.  The  prediction  Interval  (or  scan  Interval) 
begins  at  tj  so  that  (tj-t^)  represents  the  history  Interval  vihlch 
contains  m^  bearings.  The  number  of  bearlngjs  from  tj  to  is 

(m  -m.)  =  p  and  the  number  from  t^  to  t^  is  N.  For  further  slmpll- 

'  J+1  j'  o  ^ 

flcatlon  without  loss  of  generality,  B(t^)  and  t^  are  both  set  to 
zero.  For  line  1, 


B  =  b^t 


(C-5-1) 


and  for  line  2, 


B  =  a2+b2t  =  B( t2)+b2( t-t^)  =  (b^-b2) t2+b2t .  (C-5-2) 


Thus , 


and 


B*ij  =  b,t,  for  tj  <  t,  «  tj  ^  , 


B 


IJ 


for  t  .  <  t .  ^  t^ 
J  1  ^  z 


(bi-b2)tz+b2ti  for  ^  t 


J+1' 


(C-5-3) 


(C-5-4) 
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Using  equation  (C-1-3) 

N  "‘j+l 

l=^j+l  1-N 

In  which  the  first  sum  Is  0  by  virtue  of  (c-5-3)  and  (c-5-4).  Intro¬ 
ducing  the  expressions  for  and 


(C-5-5) 

f 

Assuming  equal  bearing  sampling  Intervals, 

t^  -  IT  and  =  NT 

In  which  T  Is  the  time  between  samples.  Ihen, 

m .  , 

J+1 

■■j  -  I  [(bj-bjjd-N)]  . 

1=N 

Now  let  1  «  K+N,  then 


'’j  “  I  [(^-^2)^] 

K=>0 
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or 


(bi-bg 

- 5 — 


) 


(C-5-6) 


Since  N  by  hypothesis,  only  (b^-bg)  can  be  negative.  Also, 

from  equations  (C-5-1,  2) 


(bi-bg)!  -  |ab 


so  that  (C-5-6)  reduces  to 


or 


(C-5-7) 


for  tj<  Vr 

Consider  now  tj  as  Indicated  In  Pig.  C-5-2.  'fhe  least-square  fit 

of  the  bearings  from  t^  to  tj  Is  represented  by  line  J  which  will  In¬ 
tersect  line  1  at  t^/S  and  line  2  at  {t.-tj^)/2.  The  equation  of  line 
J  is 


t 


|(bi-b2)tz+b2tjJ+ ^  {(bi-bg)  }•  (C-5-8) 

^  J  J 


Using  this  equation  for  B  and  the  equation  of  line  2  for  B.  .,  equa¬ 
tion  (C-1-3)  reduces  to 


■| 

.1 

'I 

I 


.  1 
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1  AB  ‘'z 


‘J  S’ 


(C-5-9) 


for  t j .  Tills  relation  reduces  to  equation  (C-5-7)  vfhen  t^  ■  t j , 

as  Indeed  It  should. 


It  will  subsequently  be  seen  that  for  small  Ab  the  ratio  ^j/°j  de¬ 
creases  with  Increasing  This  merely  reflects  the  results 

of  extrapolation  and  should  be  expected.  Suppose  a  zlg  has  occurred 
In  the  Interval  that  P^j  was  not  of  sufficient  magni¬ 

tude  to  conclude  that  a  zlg  occurred.  Extrapolation  Is  terminated  at 
tj^l  whence  all  of  the  bearing  data  prior  to  Is  least-square 

fitted,  and  the  Interval  (tj^2“^j+l^  scanned  using  this  latter 
least-square  fit,  and  so  on,  for  the  next  Interval.  We  would  like  to 
determine  the  behavior  of  Tj  In  these  subsequent  scanning  Intervals. 
The  situation  Is  Illustrated  In  Pig.  C-5-3.  given  by  equation 

(C-5-9).  Thus, 


Jj+K  _  ^J+K-1^^1+K+l‘'‘^J+K^^^J+K+l“'‘^z'^'^^ 


'j+K-1  ^J+K(VK-tj+K-l)(tj+K-'^z-^T) 


If  t 


J+1+1 


tj^j^  «•  At  for  all  1, 


J+K-1 


which  for  T  small  compared  to  reduces  to 


Vj+K“^z/ 


J+K-l 


■J+K  Vj+K 


(C-5-10) 
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Fig.  C-5-3  Behavior  of  Linear  P.  In  Succeeding 
Scan  Intervals  ^ 
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Since  for  K  >  1, 


J+K 


^J+K-1 


1 

for 

>  At 

-  j 

'< 

1 

1 

for 

^z 

<  At 

1 

1  " 

1 

for 

^z 

-  At 

(C-5-11) 


It  can  thus  be  concluded  that  increases  when  we  progress  to  succeed¬ 
ing  scan  intervals  whenever  >  At,  and  conversely. 


The  value  of  any  rapidly  approaches  its  limiting  value  for  a 

fixed  scan  Interval  as  gets  large.  The  limiting  value  is 


limit 
K  — ^  00 


J+K 


1  |A  B  I 
?  T 


t^At. 


(C-5-12) 


Approximate  expressions  for  AB  for  long  ranges  can  easily  be  obtained 
for  specific  geometries.  We  assume  range  to  be  approximately  constant. 
For  Pig.  C-5-^ 


AB  «=  -(V^^^slnC)A. 


and  for  Pig.  C-5-5 


AB  «  V^^^(l  -  coaOA- 


For  a  speed  zlg. 


AB  «=  AV 


(2) 


/R. 


(C-5-13) 


(C-5-14) 


(C-5-15) 


These  equations  and  the  equations  for  T.  were  used  to  obtain  Pies 
C-5-6  and  C-5-8.  ^  ’ 
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Fig.  C-5-4  ZIg  Geometry  No.  1  ^ 

3 

3 


OWN -SHIP  TRACK 


Fig.  C-5-5  ZIg  Geometry  No.  2 
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The  quadratic  bearing  extrapolation  expressed  by  equation  (C-l-l)  has 
the  disadvantage  of  having  a  much  larger  Oj  associated  with  It.  Addi¬ 
tional  difficulties  arise  when  the  zlg  is  not  detected  In  the  Interval 
In  xvhlch  It  occurred.  This  difficulty  arises  from  the  tendency  of  the 
quadratic  least-square  fit  to  follow  the  bearlng-tlme  curve  "knee" 
caused  by  the  zlg.  The  situation  Is  qualitatively  Illustrated  in  Pig. 
C_5_9  where  the  parabolic  curve  represents  a  least-square  fit  of  the 

bearing  data  from  t  to  t  . , , .  It  appears  that  linear  extrapolation 

o  j+± 

Is  superior  to  quadratic  extrapolation  for  two  reasons:  1)  for  t^ 
near  tj,  Tj  for  long  ranges  Is  approximately  the  same  for  both,  but 
the  Oj  associated  with  linear  extrapolation  Is  considerably  less;  and 
2)  zlgs  not  detected  In  the  Interval  In  which  they  occurred  have  a 
greater  probability  of  being  detected  In  subsequent  scans  with  linear 
extrapolation.  If,  however,  the  Interval  of  validity  of  the  linear 
equation  Is  comparable  to  zlg  execution  time,  as  might  occur  at  short 
ranges,  the  quadratic  method  might  be  better  than  the  linear  method. 

A  detailed  study  of  this  situation  is  needed. 


6.  Theoretical  Zlg  Probabilities 

Graphs  of  Tj  and  the  associated  Oj  based  on  the  approximate  equations 
of  the  previous  sections  are  presented  in  this  section.  Only  a  few 
geometries  are  considered.  Zlg  probability  tables  accompany  each 
graph;  the  tables  correspond  to  the  graphs  bearing  the  same  numbers. 


Figs.  C-6-1,  C-6-2,  and  the  tables  Illustrate  the  conclusion  arrived 
at  In  the  previous  section.  The  zlg  was  assiuned  to  occur  at  the  be¬ 
ginning  of  the  scan  Interval. 

Figs.  C-6-3  through  0-6-5  and  their  tables  Illustrate  zlgs  which  have 
small  effect  on  bearing-rate.  It  Is  very  unlikely  that  these  zlgs 
would  be  detected  using  the  quadratic  method. 
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In  all  of  the  cases  considered,  speed  reductions  during  the  zigs  would 
tend  to  increase  r^.and  consequently  advance  detection  time.  Finite 
zig  times  would,  of  course,  retard  detection  time. 
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n  LINEAR  VS.  TIME 


SCAN  TIME  IN  SECONDS 

Fig.  C-6-2  Ij  Linear  Versus  Time 
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Fig.  C-6-3  r  Linear  Versus  Time  for  Two- 
Mlnute  Scons  -  20  Nautical  Mile 


2IG  TIME  IN  MINUTES 


Fig.  C-6-4  r  Linear  Versus  Time  for  Two- 
Minute  Scons  -  15  Nautical  Mile 
Target 


ZIG  TIME  IN  MINUTES 


CONFID 


Fig.  C-6-5  r  Linear  Versus  Time  for  Two- 
Mlnute  Scons  -  10  Nautical  Mile 
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